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ABSTRACT 


This  dissertation  examines  the  electronic  structure  and  thermochemistry  of  low- 
valent  aluminum  clusters  that  may  serve  as  precursors  for  new  energetic  materials. 
Clusters  such  as  AlsoCp  n  (Cp*=CsMe5)  have  theoretical  heats  of  combustion  more  than 
twice  that  of  common  high  explosives  and  potentially  faster  combustion  kinetics  than  bulk 
metals.  The  tetrameric  aluminum  cluster  ALtCp*4  is  a  prototypical  monovalent  aluminum 
compound,  and  a  potential  precursor  for  these  larger  metalloid  clusters.  The  synthesis  of 
AI4R4  (R=C5Me4Pr,  CsM^iPr),  two  clusters  similar  to  Al4Cp  4,  was  recently  reported  and 
the  effect  of  their  increased  steric  bulk  is  discussed  here.  Experimental  results  and  density 
functional  theory  (DFT)  analysis  show  that  these  clusters  are  enthalpically  more  stable 
than  the  Cp  variant,  due  primarily  to  non-covalent  interactions  (NCIs)  across  ligand 
groups.  These  NCIs  show  how  ligand  steric  bulk  can  add  stability  to  tetrameric  clusters  in 
addition  to  low-valent  metal  bonding.  Similar  calculations  are  performed  on  seven  other 
homoleptic  Cp-related  clusters  with  varying  levels  of  steric  bulk.  DFT  results  are  used  to 
predict  monomer/tetramer  equilibrium  for  all  clusters  and  show  trends  counter  to 
expectations  regarding  the  role  of  ligand  bulk.  This  equilibrium  could  be  an  important 
component  in  determining  the  viability  of  a  cluster  as  a  precursor  for  larger  clusters. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  addition  of  aluminum  particles  to  high  explosives  is  a  common  method  of 
improving  overall  energy  production  and  can  enhance  blast  effects  [1].  PBXN-109,  for 
example,  is  a  common  general  purpose  high  explosive  that  is  composed  of  a  mix  of 
hexogen  (RDX),  aluminum,  binder  and  plasticizer.  The  addition  of  aluminum  to  an 
explosive  formulation  generally  increases  the  overall  blast  impulse,  but  reduces 
detonation  velocity  and  pressure  [2].  To  understand  the  contribution  aluminum  makes  to 
the  improvement  of  high  explosives,  a  brief  overview  of  detonation  is  necessary. 

Detonation  is  a  process  involving  a  supersonic  pressure  wave  perpetuated  by  an 
ongoing  chemical  reaction.  For  example,  the  decomposition  of  RDX  into  N2,  H2O,  CO 
and  CO2  releases  thousands  of  kilojoules  of  energy  per  mole,  creating  heat  and  pressure 
that  will  perpetuate  the  chemical  reaction  and  can  drive  metal,  a  blast  wave  or  otherwise 
do  work  external  to  the  explosive.  The  detonation  wave  moves  several  kilometers  per 
second  and  is  driven  by  a  reaction  zone  that  is  on  the  order  of  a  few  millimeters  or  less  in 
width  [3].  Only  energy  released  inside  the  reaction  zone  will  contribute  to  detonation 
wave  propagation  and  its  peak  pressure.  This  requires  the  reaction  time  of  the  high 
explosive  to  be  less  than  a  microsecond — otherwise  the  detonation  wave  and  reaction 
zone  will  pass  over  the  molecule  before  it  has  had  time  to  release  energy.  Aluminum 
burns  too  slowly  to  contribute  to  the  detonation  wave,  though  depending  on  the 
confinement  of  the  charge,  post-detonation  combustion  of  aluminum  can  still  contribute 
to  brisance  [4],  The  ignition  time  for  aluminum  varies  depending  on  the  size  and  type  of 
particle,  but  is  on  the  order  of  just  under  a  millisecond  to  tens  of  milliseconds  for  typical 
micrometer  to  nanometer  scale  particles  included  in  explosive  fonnulations  [5],  [6].  By 
the  time  the  aluminum  particles  have  ignited,  the  reaction  zone  has  passed.  For  this 
reason,  aluminum  is  generally  treated  as  an  inert  material  from  a  detonation  perspective. 
Figure  1  shows  reaction  times  for  RDX  and  a  typical  micrometer-sized  aluminum  powder 
compared  to  the  reaction  times  required  to  influence  key  capabilities  of  a  high  explosive. 
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Figure  1 .  Reaction  times  for  RDX  and  aluminum  powder  and  influence  times 

for  explosive  effects.  Adapted  from  [7]. 


Even  though  aluminum  does  not  help  sustain  a  detonation  wave,  it  is  still 
commonly  added  to  the  current  stock  of  military-grade  high  explosives  like  RDX, 
octogen  (FIMX)  and  hexanitrohexaazaisowurtzitane  (CL-20).  The  aluminum  typically 
contributions  to  blast  effects  and  has  a  relatively  high  energy  density.  Despite  more  than 
a  century  of  development,  the  energy  density  of  high  explosives  has  improved  slowly. 
Figure  2  shows  a  rough  progression  from  trinitrotoluene  (TNT)  to  CL-20,  the  most 
powerful  commodity  high  explosive  in  wide  scale  production  [8].  TNT  equivalence  for 
the  organic  compounds  was  detennined  by  the  chemical  equilibrium  code  Cheetah  and 
compares  total  energy  released  by  detonation.  Also  shown  is  bulk  aluminum.  For 
aluminum,  the  total  energy  release  by  detonation  for  TNT  is  compared  to  the  heat  of 
combustion  of  aluminum  metal  oxidized  to  form  aluminum  oxide  (AFO3),  83  kJ/cm’  [9]. 
This  represents  a  theoretical  maximum  energy  release  for  aluminum  with  infinitely 
available  oxidizer,  and  shows  the  tremendous  potential  of  metalized  explosives. 
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Explosives  are  listed  based  on  the  year  they  were 
first  synthesized. 

Figure  2.  TNT  equivalence  of  organic  explosives  and  bulk  aluminum. 

Reducing  the  timescale  of  aluminum  reaction  would  represent  a  significant 
advance  in  high  explosive  development,  allowing  the  high  energy  density  of  the  metal  to 
be  incorporated  into  detonation  and  metal-throwing  effects.  While  new  organic 
explosives  (e.g.,  octanitrocubane  [10]  and  di-l,2,3,4-tetrazine  tetraoxide  (DTTO)  [11], 

[12] )  still  hold  promise  for  significant  improvements  in  organic  explosive  performance, 
their  total  energy  density  is  unlikely  to  approach  that  of  heavily  aluminized  formulations 

[13] .  If  the  kinetics  of  the  aluminum  combustion  can  be  rapid  enough  to  react  in  or  close 
to  the  reaction  zone,  significant  enhancements  beyond  organic  explosives  could  be 
realized. 

1.  Combustion  of  Aluminum 

Aluminum  is  a  group  13  metal  with  a  natural  oxidation  state  of  +3.  It  exists 
naturally  in  minerals  in  the  +3  oxidation  state  and  the  pure  metal  form  is  only  obtained 
through  industrial  processes.  Bulk  aluminum  readily  reacts  with  oxygen  forming  the 
ceramic  AFCfi  releasing  about  31  kJ/g  of  energy  [14],  By  comparison,  1  gram  of  CL-20 
releases  about  8  kJ/g  of  energy  [15].  Because  of  the  increased  energy  density  of 
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aluminum  compared  to  organic  explosives,  it  is  routinely  used  as  an  energetic  additive  in 
many  explosive  compounds.  Typically,  aluminum  is  added  in  the  form  of  small 
manufactured  particles  ranging  in  size  from  about  100  nanometers  to  a  few  hundred 
microns.  Their  relatively  long  ignition  time  and  burn  time,  however,  prevent  the  energy 
released  during  oxidation  from  participating  in  the  detonation  reaction. 

A  typical  aluminum  micro-particle  is  not  purely  bulk  aluminum.  Though  usually 
synthesized  in  an  argon  atmosphere,  particles  are  eventually  exposed  to  air.  Once 
exposed,  the  surface  of  a  particle  will  immediately  oxidize  and  produce  a  thin  ceramic 
shell  around  the  interior  aluminum  known  as  a  passivation  layer.  The  thickness  of  the 
passivation  layer  varies,  but  is  usually  around  a  few  nanometers.  Figure  3  shows  a 
transmission  electron  microscope  (TEM)  image  for  a  roughly  100-nanometer  diameter 
aluminum  particle.  The  passivation  layer  represents  the  first  barrier  to  efficient 
combustion  of  an  aluminum  particle.  Before  the  bulk  aluminum  core  of  a  particle  can  be 
oxidized,  the  oxidizer  must  penetrate  the  aluminum  oxide  shell.  This  is  generally 
accomplished  in  one  of  two  ways:  either  by  the  diffusion  of  the  oxidizer  and  aluminum 
metal  through  the  heated  aluminum  oxide  shell,  or  due  to  the  fracturing  of  the  aluminum 
oxide  shell  as  a  result  of  the  expansion  of  the  expanding,  molten  aluminum  core. 
Regardless  of  the  mechanism,  delay  in  the  ignition  of  the  aluminum  particle  will  result. 
This  delay  depends  on  a  number  of  factors,  but  at  its  shortest  is  just  under  a  millisecond 
[6].  This  delay  is  sufficient  to  prevent  aluminum  oxidation  from  participating  in  the 
detonation  process.  Efforts  to  passivate  aluminum  with  something  other  than  aluminum 
oxide  are  ongoing  and  success  has  been  achieved  in  the  area  of  self-assembled 
monolayers  [16].  Jouet  et  al.  successfully  produced  organic-coated,  oxide-free  aluminum 
nanoparticles  that  were  air  stable  [16].  While  this  approach  potentially  removes  a  barrier 
to  ignition  of  aluminum  particles,  it  reduces  the  overall  aluminum  content  compared  to 
traditionally  passivated  aluminum  particles.  This  is  because  the  molar  mass  of  the 
passivation  agents  used  are  25-30  times  greater  than  that  of  bulk  aluminum  [16]. 
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Figure  3.  TEM  image  of  aluminum  particle  showing  aluminum  oxide 

passivation  layer.  Source:  [17]. 

While  production  of  oxide-free  aluminum  particles  would  reduce  ignition  times 
and  potentially  provide  an  oxidizing  agent  for  the  aluminum,  there  is  a  second  barrier  to 
efficient  burning  of  bulk  aluminum.  The  boiling  point  of  aluminum  oxide  is  about  500  K 
higher  than  that  of  aluminum.  During  particle  combustion,  gas-phase  aluminum  and 
oxygen  gas  create  aluminum  oxide.  The  aluminum  oxide  tends  to  condense  and  form  a 
cap  which  attaches  to  the  remaining  aluminum  in  the  particle  [18].  This  oxide  cap  further 
complicates  and  inhibits  the  combustion  process  and  ultimately  slows  the  burn  rate  of  the 
particle.  Figure  4  shows  a  70-micrometer  diameter  aluminum  particle  about  three 
milliseconds  after  ignition.  The  smooth  structure  on  the  right  is  aluminum  oxide.  This 
oxide  cap  prevents  aluminum  particles  from  combusting  in  a  manner  similar  to 
hydrocarbon  droplets,  which  have  a  burn  time  proportional  to  the  diameter  of  the  droplet 
squared.  Efforts  to  model  aluminum  particle  combustion  are  ongoing,  but  what  is  relevant 
to  this  work  is  that  bum  time  is  generally  on  the  order  of  milliseconds.  While  burn  times 
from  aluminum  microparticles  have  been  measured  in  the  hundreds  of  milliseconds, 
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Beckstead  points  out  that  the  vast  majority  of  the  available  energy  is  released  within 
about  the  first  five  milliseconds  [19].  Work  by  Dreizin  corroborates  this,  showing  that  for 
90  micrometer  aluminum  particles,  burn  times  at  3,000  K  were  about  5-10  milliseconds 
[20].  Reducing  particle  size  to  the  nanometer  regime  does  improve  bum  times.  Shock 
ignition  of  100  nanometer  aluminum  particles  had  burn  times  of  just  under  600 
microseconds  [6].  This  is  almost  ten  times  faster  than  a  typical  micron-sized  particle,  but 
still  significantly  slower  than  the  reaction  time  of  a  molecule  like  RDX,  less  than  one 
microsecond.  There  is  a  drawback  to  using  small,  nanometer-sized  particles,  however.  As 
size  decreases,  the  aluminum  oxide  passivation  layer  makes  up  a  larger  and  larger 
percentage  of  the  total  mass  of  the  particle.  The  impact  varies  depending  on  the 
manufacturing  process  and  the  age  of  the  particle,  but  in  general,  particles  less  than  about 
50  nanometers  are  composed  of  no  more  than  about  40%  of  aluminum  by  mass  [16]. 


Figure  4.  Aluminum  particle  quenched  2.5-3  milliseconds  after  ignition. 

Source:  [18]. 


Typical  aluminum  particles  used  in  explosive  compounds  ignite  and  bum  long 
after  the  detonation  wave  has  passed.  Even  advances  in  aluminum  particle  production, 
like  size  reduction  and  reduction  or  even  elimination  of  the  passivation  layer,  still  put 
energy  release  2-3  orders  of  magnitude  too  slow  to  contribute  to  detonation  wave 
propagation.  Additionally,  these  ignition  and  bum  time  calculations  assume  a  readily 
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available  supply  of  oxidizer,  typically  O2,  in  the  gas  phase.  This  is  not  likely  to  be  the 
case  when  the  particles  are  packed  into  an  explosive  compound.  This  work  explores  an 
alternative  avenue  to  harnessing  the  full  energetic  potential  of  bulk  aluminum.  Ongoing 
research  into  small  organometallic  aluminum  clusters  has  shown  that  clusters  with  large 
numbers  of  low  valent  aluminum  atoms  can  be  synthesized.  While  considerable  synthesis 
challenges  remain,  their  relatively  small  size  and  potential  to  include  oxidizer  within  the 
surrounding  ligands  could  allow  them  to  oxidize  on  time  scales  similar  to  those  of 
traditional  organic  explosives. 

2.  Low-Valent  Aluminum  Clusters 

Monovalent  aluminum  clusters  have  received  considerable  attention  in  recent 
years,  following  the  initial  synthesis  of  the  monovalent  aluminum  cluster  AfiCp  4  (Cp  = 
C'sMcs)  [21].  Investigations  on  this  and  similar  low  valence  aluminum  clusters  have 
continued;  AfiCp*4  has  served  as  an  aluminum  atom  source  in  the  fonnation  of 
nanoparticles,  and  AlCp  units  have  proven  to  be  an  efficient  ligand  in  many  transition 
metal  complexes  [22]— [25].  Larger  clusters  with  a  significant  number  of  low-valence 
aluminum  atoms  have  also  been  investigated  [26]-[30].  These  systems  are  of  interest  for 
a  variety  of  applications,  including  use  as  energetic  material  [31].  The  unique  structure  of 
these  clusters,  with  aluminum  atoms  in  the  central  core  surrounded  by  organic  ligands, 
may  allow  for  fast  metal  combustion  kinetics  while  avoiding  the  rapid  sintering  and 
native  oxide  layer  that  limit  traditional  aluminum  nanoparticle  combustion  [32].  The 
close  proximity  of  exterior  ligands  to  the  low-valent  metal  core  may  offer  opportunities 
for  very  rapid  mass  transfer  between  the  metal  and  oxidizing  functional  groups  on  the 
ligand.  The  energetic  potential  of  Al4Cp  4  and  a  similar,  larger  cluster,  Al50Cp  12,  were 
shown  to  potentially  have  volumetric  heats  of  combustion  approaching  60%  of  bulk 
aluminum  [31].  This  is  especially  impressive  considering  the  density  of  these  clusters  is 
less  than  half  of  that  of  bulk  aluminum  [31].  Figure  5  shows  a  comparison  of  calculated 
heats  of  combustion  by  mass  of  both  compounds  compared  to  RDX  and  a  generic 
aluminized  plastic  bonded  explosive,  PBX. 
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Figure  5.  Fleat  of  combustion  by  mass  of  two  organometallic  clusters  and  two 

typical  organic  explosives. 

Compared  to  aluminum  particles  currently  used  in  explosives,  these  low  valent 
aluminum  clusters  are  expected  to  bum  significantly  faster  than  traditional  bulk  metals. 
There  is  no  oxide  layer  surrounding  these  particles  and  an  oxidizer  could  potentially  be 
integrated  into  the  surrounding  organic  ligand.  All  that  is  between  the  aluminum  core  and 
the  oxidizing  agent  is  a  stabilizing  organic  ligand,  which  is  only  a  few  bond  lengths  thick. 
Additionally,  these  clusters  are  significantly  smaller  than  the  smallest  aluminum  particles 
used  in  explosives.  Figure  6  shows  an  aluminum  nanoparticle  with  a  diameter  of 
25  nanometers  and  an  oxide  layer  three  nanometers  thick.  The  oxide  layer  represents 
about  56%  of  the  total  volume  of  the  particle.  To  the  right  are  to-scale  depictions  of 
Al5oCp*i2  and  RDX.  This  visual  demonstrates  not  only  the  increased  impact  of  the  oxide 
layer  on  small  aluminum  nanoparticles,  but  also  the  dramatic  size  difference  between 
these  particles  and  even  large  monovalent  aluminum  clusters. 
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Figure  6.  A  to-scale  depiction  of  a  25-nanometer  diameter  aluminum  particle, 

Al50Cp*i2  and  RDX. 

Currently  there  are  seven  synthesized  and  structurally  characterized  tetrameric 
clusters  containing  a  tetrahedral,  monovalent  aluminum  core  [33]— [39] .  Of  these  seven, 
only  two  are  homoleptic  complexes  stabilized  by  a  variant  of  the  cyclopentadienyl  ligand: 
Al4Cp*4  and  Al4Cp4  (Cp  =C5Me4H).  The  tetramer  of  the  unsubstituted  variant,  Al4Cp4 
(Cp=CsH5),  has  been  observed  in  solution  but  never  isolated  [38].  Work  on  these  clusters, 
as  well  as  other  low  valence  aluminum  clusters,  is  ongoing;  as  the  synthesis  and  isolation 
of  these  materials  is  difficult  and  the  final  products  are  often  air-sensitive,  density 
functional  theory  (DFT)  analysis  has  been  very  helpful  in  understanding  the  behavior  of 
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these  clusters  [40]— [46].  These  tetrameric  clusters  are  of  importance  because  they 
represent  potential  precursors  for  larger  metalloid  clusters  [47].  Of  chief  concern  is  the 
monomer  to  tetramer  ratio  of  these  clusters  in  equilibrium.  A  high  ratio  of  monomers  to 
tetramers  is  desired  for  a  large  cluster  synthesis  approach  pursued  by  Curtis  Johnson  and 
David  Morris  at  Naval  Air  Warfare  Center  (NAWC),  China  Lake.  The  approach  takes 
advantage  of  two  known  reactions,  shown  in  the  top  two  positions  in  Figure  7.  The 
bottom  reaction  in  the  Figure  is  the  approach  of  Johnson  and  Morris  [48].  The  first 
reaction  demonstrates  the  ability  of  the  AlCp  variant  to  react  with  the  trivalent  species, 
AlCfiH.  The  second  reaction  demonstrates  the  ability  of  potassium  graphite  to  reduce 
trivalent  AlCfiH  to  bulk  aluminum  metal.  By  using  a  low  molar  amount  of  potassium 
graphite,  to  prevent  reduction  of  AlCp  monomers  and  tetramers,  the  hope  is  that  bulk 
metal  produced  by  the  reduction  of  A1CLH  in  combination  with  available  AlCp 
monomers  in  solution  will  lead  to  large,  metalloid  clusters  analogous  to  AlsoCp*^.  A  high 
ratio  of  monomers  to  tetramers  for  Al4Cp  4  and  analogous  compounds  would  provide 
more  potential  building  blocks  for  larger  clusters. 


[AlCp*]4  +  A1CLH-2THF 

dark  grey  powder  +  unidentified  product  +pale  yellow  solution 

Toluene 

A1CLH-2THF 

(0.1-2)  KC8 

Toluene 

A1  metal 

[AlCp*]4  +  AlCl2H-2THF 

0.1  KCg 

large  metalloid  clusters 

Toluene 

Figure  7.  Large  cluster  synthesis  approach. 


Until  recently,  monomers  in  solution  were  only  detectable  for  the  Cp*  variant  with 
a  ratio  of  about  one  monomer  for  every  500  tetramers  at  room  temperature.  Two  new 
aluminum  clusters  with  substituted  cyclopentadienyl  ligands  have  been  recently 
synthesized  by  Dennis  Mayo  and  Rebecca  Wilson  at  the  Naval  Surface  Warfare  Center 
(NSWC),  Indian  Head  Explosive  Ordnance  Disposal  Technical  Division  (EODTD),  both 

with  detectable  monomers  in  solution.  The  properties  of  these  new  systems  are 
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considered  in  great  detail  in  following  sections.  Both  clusters  had  monomer  to  tetramer 
ratios  greater  than  the  Cp  variant.  While  these  systems  failed  to  crystallize  and 
unambiguous  structural  data  could  not  be  obtained  experimentally,  a  combination  of 
variable-temperature  27A1  nuclear  magnetic  resonance  (NMR)  experiments  and  DFT 
calculations  described  later  strongly  indicate  these  are  tetrameric  Al4Cp*Pr4 
(Cp  P -(AMe4Pr)  and  Al4Cp*lPr4  (Cp*lPr=C5Me4iPr)  and  their  respective  monomers.  These 
additional  clusters  provide  key  experimental  thermochemical  data  that  will  ultimately 
help  in  developing  a  computational  means  of  predicting  monomer  to  tetramer  ratio  for 
theoretical  tetrameric  clusters.  Computational  approaches  to  electronic  structure  and 
thermochemistry  analysis  are  varied,  however,  and  different  methodologies  can  lead  to 
substantially  different  results.  Before  discussion  of  computational  approaches  to 
analyzing  low-valence  aluminum  metalloid  clusters,  it  is  important  to  understand  the 
origins  and  structure  of  modem  DFT. 


B.  DENSITY  FUNCTIONAL  THEORY 


DFT  provides  an  avenue  for  determining  the  ground  state  of  a  many  electron 
system.  In  theory,  exact  solutions  can  be  obtained,  however  in  practice  this  is  not  the 
case.  DFT  has  its  roots  in  the  Schrodinger  equation  which  is  analytically  solvable  for 
single  electron  systems,  but  becomes  extremely  complex  for  systems  with  even  just  a  few 
electrons.  For  a  system  with  n  electrons  and  M  nuclei  the  many-body  Hamiltonian  is: 


HB0  —  TP  +  VpN  +  VPP  +  V, 


eN 


vee  1  VNN 
n  M 


-£>ZZ- 

i= 1  i= 1  1=1 


n  n 

V*1  V’1 


R, 


+ 


Z_t  \r  —  r  | 


+  const. 


i= 1  j>i 


(1) 


Here  we  use  atomic  units  where  the  mass  of  the  electron,  its  charge,  h  and  47T£0 
are  all  set  to  unity.  The  atomic  unit  of  energy  is  called  the  Hartree  and  is  -27.21  eV.  This 
equation  also  makes  use  of  the  Born-Oppenheimer  (BO)  approximation  which  assumes 
the  nuclei  are  motionless.  This  eliminates  the  kinetic  energy  tenn  for  the  nuclei  and  will 
reduce  the  nuclei-nuclei  potential  to  a  constant.  To  detennine  the  wave  functions  for  the 
ground  state  of  this  system  we  need  to  solve  the  time-independent  Schrodinger  equation: 
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HBOxi,(x1,...,xn)  =  E0'¥(x1,...,xn).  (2) 

The  vectors  xL  denote  the  coordinates  of  the  /th  electron  (rL)  as  well  as  the  spin  of 
the  7th  electron  (up  or  down).  This  is  a  second  order  differential  equation  with  4N 
variables.  Exact  solutions  are  beyond  reach.  A  key  property  of  the  Hamiltonian,  however, 
is  that  when  it  acts  on  the  true  wavefunction  the  lowest  possible  energy  is  produced.  This 
allows  the  use  of  the  variation  principle  where  an  initial  guess  of  the  wavefunction  is 
continuously  varied  until  the  energy  is  minimized. 

1.  Hartree-Fock  Theory 

Hartree-Fock  (HF)  theory  is  an  early  attempt  to  arrive  at  the  wavefunction  of  a 
many-body  system  developed  in  the  1920s  and  1930s.  HF  theory  makes  use  of  a  very 
simple  guess  for  the  structure  of  a  many  body  wavefunction  and  then  makes  use  of  the 
variation  principle  to  find  an  approximation  of  the  actual  wavefunction.  First,  the  theory 
assumes  that  all  electrons  in  the  system  are  non-interacting.  The  wavefunction  would 
then  just  be  a  combination  of  single-electron  wavefunctions.  Simply  multiplying  single 
electron  wavefunctions  together  would  not  produce  the  anti-symmetry  required  for 
wavefunctions  of  fermions,  so  the  appropriate  form  is  a  Slater  determinant  (SD): 

1  i/fiOi)  -  i) 

^SnO'T'  ■■■  ’  )  ^[tl\  "  ''  •  (3) 

'  Ipifrn)  ...  lpn(xn) 

The  single  electron  wavefunctions,  xpi,  are  called  orbitals  and  make  up  an  orthonormal 

set.  When  the  three  non-constant  terms  of  HB0,  Eq.  (1),  act  on  fSD,  four  terms  are 

produced.  The  first  two  can  readily  be  identified  as  the  kinetic  energy  and  the  electron- 

nucleus  interaction  terms.  The  third  and  fourth  terms  are  both  associated  with  the  two- 

particle  summation  that  represents  electron-electron  interactions,  Vee .  They  are  known  as 

the  Coulomb  (/)  and  exchange  ( K )  terms  respectively.  The  reason  both  exist  is  purely  a 

consequence  of  the  antisymmetric  nature  of  the  Slater  determinant  [49].  J  represents  the 

self-interaction  of  the  electron  density  cloud  and  the  potential  is  considered  local  [49], 

[50].  K  represents  self-repulsion  and  the  potential  is  considered  non-local  [49],  [50]. 

Because  all  terms  contain  only  combinations  of  single  particle  wavefunctions,  they  can 

be  readily  evaluated  and  the  ground  state  energy  calculated.  The  orbitals  are  then  varied, 
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preserving  their  orthononnality,  until  the  energy  is  minimized.  Despite  the  simple  initial 
guess  at  the  true  many  particle  wavefunction,  HF  theory  can  produce  ground  state 
energies  that  are  as  much  as  99.9%  accurate  [50].  The  energy  the  HF  theory  fails  to 
capture  is  called  correlation  energy,  Ec,  and  can  be  quite  large  for  even  relatively  small 
molecules. 

2.  Density  Functional  Theory 

DFT  is  an  effort  to  improve  upon  the  success  of  HF  theory  by  finding  a  route  to 
the  ground  state  energy  with  only  knowledge  of  the  electron  density,  p(r).  In  theory,  all 
we  need  to  fully  determine  HB0  is  knowledge  of  the  number  of  electrons  and  the  potential 
they  are  subject  to.  This  potential  is  defined  by  the  static  positions  of  the  atomic  nuclei 
and  the  number  of  electrons  can  be  determined  by  integration  of  p(r )  over  all  space. 
Figure  8  shows  a  plot  of  p(r )  for  the  Cp"  anion.  The  positions  of  the  atomic  nuclei  can 
clearly  be  seen  as  cusps  of  p (r) .  The  cusps  are  related  to  the  nuclear  charge  of  the 
nucleus  by  the  equation: 


where  p(r)  is  the  spherical  average  of  p(r)  [49].  The  electron  density  allows  us  to 
determine  not  only  nuclei  position,  but  also  atomic  charge. 
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Figure  8.  Electron  density  plot  of  Cp"  anion. 


The  final  key  element  is  a  theorem  by  Hohenberg  and  Kohn  which  states  that  any 
ground  state  configuration  of  electrons  is  uniquely  determined  by  the  external  potential 
[51].  Stated  another  way,  a  certain  configuration  of  nuclei  surrounded  by  a  certain 
number  of  electrons  will  give  rise  to  one  and  only  one  p  (r) .  That  means  there  is  a  unique 
ground  state  energy  and  wavefuntion  associated  with  HB0,  which  we  can  get  from  p(r). 
This  is  because  p(r )  can  provide  nuclear  positions  and  charges,  and  integration  of  p(r) 
can  provide  the  number  of  electrons.  All  the  information  we  need  to  determine  HB0,  and 
therefore  the  ground  state  energy,  must  be  available  within  the  electron  density  and  so 
there  must  be  some  energy  functional  E\p(r)\  which  produces  the  correct  ground  state 
energy  if  given  the  correct  electron  density.  The  problem  is  reduced  from  a  4N 
dimensional  second  order  differential  equation  to  a  functional  with  three  special 
variables. 

The  form  of  this  energy  functional  is  not  clear.  The  electron-nucleus  interaction 
potential,  VeN,  can  be  readily  calculated  from  p(r),  but  the  kinetic  energy  and  the 
electron-electron  interaction  potential  cannot  [50].  Ultimately  the  form  of  the  energy 
functional  is  motivated  from  FIF  theory.  About  80-90%  of  the  electron-electron  potential 
is  captured  in  the  FIF  term,  /,  which  can  be  calculated  from  p(r)  [50].  The  kinetic  energy 
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of  non-interacting  particles,  Ts,  is  indirectly  a  functional  of  p(r)  and  is  a  large  portion  of 
the  true  kinetic  energy  [50].  The  energy  functional  is  then: 

E[p]  =  Ts[p]  +  VeN[p]  +J[p]  +  Exc[p]  +  Vnn>  (5) 

where  VNN  is  the  constant  nucleus-nucleus  potential  and  the  exchange-correlation  energy, 

EXM,  is: 

ExM  =  K[p]  +  Ec[p]  +  (T[p]  -  Ts[p ]),  (6) 

where  K[p]  is  the  HF  exchange  term  and  Ec[p\  is  the  correlation  energy,  the  energy 
missing  from  HF  theory.  The  fonn  of  Exc[p]  is  unknown,  but  if  it  was  known  then  E[p] 
would  be  exact  [49]. 

The  best  modern  DFT  can  do  is  to  continually  create  better  and  better 
approximations  of  Exc  lp\.  The  two  most  common  approaches  are  based  on  a 
hypothetical  uniform  electron  gas.  The  local  density  approximation  (LDA)  uses  the 
known  form  of  the  exchange  energy  for  a  unifonn  electron  gas  and  very  accurate 
simulation  data  to  motivate  the  remaining  portion  of  Exc[p ]  [49].  The  resulting 
expression  for  Exc  [p]  is  then  applied  to  each  differential  volume  element  of  the  actual 
(non-unifonn)  electron  density  and  integrated  to  obtain  the  total  exchange-correlation 
energy  [50].  The  generalized  gradient  approximation  (GGA)  improves  upon  LDA  by 
including  terms  containing  the  gradient  of  the  actual  (non-uniform)  electron  density  [49]. 
Attempts  to  improve  upon  GGA  functionals  with  higher  derivatives  and  kinetic  energy 
terms  has  been  tried  as  well.  These  functionals  are  called  meta-GGAs  [50].  Ultimately  all 
of  these  approximations  are  theoretically  motivated,  but  include  empirical  fitting 
parameters  to  try  and  produce  the  most  accurate  results  [50].  They  are  also  considered 
local.  They  attempt  to  represent  Exc  [p],  which  contains  the  non-local  HF  exchange 
energy,  K,  with  expressions  that  only  depend  on  the  density  and  its  gradient  at  a  specific 
point  [49]. 

In  addition  to  LDA,  GGA  and  meta-GGA,  another  approach  that  is  common  in 
thermochemistry  calculations  is  the  use  of  hybrid  functionals.  Hybrid  functionals  mix 
elements  of  HF  theory  with  traditional  DFT.  The  important  element  they  include  from  HF 
theory  is  the  exact  exchange  term,  K  [49].  The  amount  varies  depending  on  the 
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formulation,  but  usually  about  20-30%  of  exact  exchange  is  included  [49].  The  inclusion 
of  HF  exchange  adds  a  non-local  element  to  the  functional  and  the  amount  is  controlled 
by  one  of  several  parameters  that  adjust  all  terms  in  the  functional.  These  parameters  are 
empirically  determined  in  order  to  produce  the  best  results  for  as  broad  a  range  of 
systems  as  possible. 

While  modem  DFT  has  enjoyed  great  success,  the  accuracy  of  functionals  is 
highly  dependent  on  the  type  of  system  investigated.  It  is  not  appropriate  to  say  that  one 
particular  functional  is  always  more  accurate  than  another,  but  in  general  GGAs  are  more 
accurate  than  LDAs  and  hybrids  are  more  accurate  than  GGAs  [49].  Of  course,  there  are 
many  types  of  functionals  within  each  group  and  in  some  cases,  for  example,  a  GGA 
functional  may  perform  better  than  a  hybrid  for  a  specific  system.  Once  a  particular 
functional  is  chosen,  then  the  ground  state  energy  of  a  specific  system  can  be 
approximated  using  the  variation  principle.  E[p ]  is  ultimately  a  functional  of  electron 
density,  but  DFT  is  based  on  the  Kohn-Sham  (KS)  approach  which  uses  orbitals  of 
individual  electrons  to  calculate  the  electron  density  [52].  These  orbitals,  called  Kohn- 
Sham  orbitals,  form  a  set,  {(p  J,  and  the  electron  density  can  be  calculated  by: 

n 

p(r)  =  ^\(Pi(r)\2. 

i 

These  KS  orbitals  are  also  used  to  directly  calculate  the  kinetic  energy  of  non¬ 
interacting  particles.  The  orbitals,  not  the  actual  electron  density,  are  varied  to  minimize 
the  energy  functional,  E[p].  They  are  constructed  from  collections  of  functions  called 
basis  sets,  discussed  in  the  next  section.  Once  the  ground  state  energy  and  electron 
density  are  known,  then  all  important  properties  of  the  system  can  be  determined.  For 
example,  the  first  derivatives  of  the  total  energy  with  respect  to  nuclei  perturbation 
provide  a  means  for  molecular  geometry  optimization  [50].  The  second  derivatives  of  the 
total  energy  can  yield  vibrational  frequencies  from  which  thermochemical  data  can  be 
calculated  via  statistical  thermodynamics  [49].  Mixed  second  derivatives  of  the  total 
energy  with  respect  to  perturbations  of  nuclear  magnetic  moments  provide  a  means  for 
determining  NMR  chemical  shifts  [49].  The  ability  of  DFT  to  detennine  these  and  many 
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other  observables  of  chemical  systems  has  made  DFT  an  essential  aspect  of  chemical 
analysis  in  general. 

3.  Basis  Sets 

The  construction  of  KS  orbitals  is  done  through  what  is  called  the  linear 
combination  of  atomic  orbitals  (LCAO).  Initially,  known  eigenfunctions  of  the  hydrogen 
atom  were  used  to  create  basis  sets  [49].  Today  there  are  four  main  approaches  to 
constructing  basis  sets,  but  the  most  common  for  treatment  of  small  molecules  and 
clusters  is  through  the  use  of  Gaussian  functions.  KS  orbitals  are  created  through  linear 
combinations  of  a  set  of  these  functions.  If  the  set  was  complete,  requiring  an  infinite 
number  of  functions,  then  the  resulting  orbitals  would  be  exact.  In  practice,  however,  the 
set  has  a  finite  number  of  functions.  The  more  functions  in  the  set,  the  more  accurately 
the  KS  orbitals  can  be  described.  This  leads  to  the  existence  of  several  different  sizes  of 
sets.  The  minimum  number  of  functions  necessary  to  describe  a  KS  orbital  is  one 
function  for  each  atomic  orbital  including  the  valence  orbitals.  Carbon,  for  example, 
would  require  five  functions,  one  for  each  of  Is,  2s,  2px,  2py  and  2pz.  A  double-zeta  basis 
set  would  have  two  functions  per  atomic  orbital;  carbon  would  require  ten  functions. 
There  are  also  triple,  quadruple  and  higher  zeta  basis  sets.  To  improve  efficiency,  split 
valence  basis  sets  were  created.  A  split  valence  double-zeta  basis  set  would  have  one 
basis  function  for  core  atomic  orbitals  and  two  functions  for  each  valence  orbital.  In  this 
case  carbon  would  require  nine  basis  functions.  It  is  also  typical  to  include  polarization 
functions  within  basis  sets.  Polarization  functions  represent  the  first  unoccupied  atomic 
orbitals  of  an  atom;  for  carbon,  these  are  the  d  orbitals.  These  functions  help  the  KS 
orbitals  distort  appropriately  in  a  molecular  environment  [49].  Just  like  DFT  functionals, 
there  are  many  types  of  basis  sets  that  are  of  different  sizes  and  constructed  in  slightly 
different  ways.  In  general,  the  use  of  a  triple-zeta  basis  set  with  polarization  is  sufficient 
to  ensure  resulting  calculations  represent  the  deficiencies  of  the  chosen  functional,  rather 
than  the  basis  set  [50].  Basis  sets  of  this  size  can  be  computationally  expensive,  however, 
so  for  large  calculations  it  may  be  necessary  to  use  smaller  basis  sets.  Smaller  basis  sets 
may  not  be  as  accurate,  but  they  can  still  produce  meaningful  results  and  capture  trends. 
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c. 


METHODOLOGY 


The  main  goal  of  this  work  is  to  understand  the  thermodynamic  properties  and 
electronic  structure  of  tetrameric,  monovalent,  aluminum  clusters  as  a  result  of  systematic 
changes  in  ligand  steric  bulk.  Of  particular  interest  is  the  ratio  of  monomers  to  tetramers 
in  solution.  Large  concentrations  of  monomers  are  likely  necessary  for  the  further 
development  of  larger  low-valent  clusters  such  as  AfioCp  12  via  the  solution-based 
synthesis  techniques  previously  discussed.  Until  recently,  experimental  thermodynamic 
data  for  monomer/tetramer  equilibrium  in  small  monovalent  aluminum  clusters  was  only 
available  for  the  Cp*  variant  via  van ’t  Hoff  analysis  of  variable-temperature  27A1  NMR 
[38].  Tetramers  of  the  Cp  and  Cp  variants  have  been  observed  in  solution,  but  both 
disproportionate  to  trivalent  aluminum  species  and  bulk  metal  before  AIR  monomers  are 
observed.  DFT  calculations  show  that  Cp1’1  and  Cp*lPr  are  bound  into  tetrameric  forms 
more  strongly  than  the  Cp*  variant,  despite  the  addition  of  steric  bulk  on  the  ligands. 

DFT  calculations  were  carried  out  to  examine  the  electronic  structure  of  these 
clusters  and  explore  changes  in  bonding  with  the  ligand  variations.  As  would  be 
expected,  the  calculated  thermochemistry  for  these  clusters  is  very  sensitive  to  the 
computational  methodology,  especially  the  choice  of  the  DFT  functional.  Previous 
calculations  gave  a  wide  range  of  results  for  tetramerization  enthalpies  and  relative 
stabilities  [31],  [38].  The  key  challenge  for  these  clusters  is  accurately  capturing  the 
dispersive  forces  present  across  ligand  groups  as  well  as  the  Al-Al  bond  behavior  in  the 
core.  The  clusters  are  combinations  of  four  separate  monomers  held  together  not  only  by 
aluminum  lone  pair  overlap  in  the  core  [31],  but  also  by  NCIs  across  neighboring  ligand 
groups.  Because  these  forces  occur  over  relatively  long  distances,  LDA,  GGA  and  meta- 
GGA  functionals,  which  are  considered  local,  are  severely  deficient  in  capturing  the 
effects  [53].  Hybrid  functionals,  which  include  the  non-local  exchange  from  HF  theory 
are  better  options,  but  in  general  they  still  cannot  reliably  describe  dispersion  [53]. 
Explicit  dispersion  corrections  can  be  calculated  by  adding  tenns  to  the  KS  energy 
functional  that  capture  the  energy  due  to  two  and  three  body  interactions.  A  common 
approach  was  developed  by  Stefan  Grimine,  called  D3  [54],  This  method  has  been  shown 
to  be  very  effective,  but  it  is  also  computationally  expensive.  Recently  work  by  Zhao  and 
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Truhlar  produced  a  suite  of  functionals  that  are  specifically  parameterized  to  handle 
compounds  with  significant  NCIs  [55].  These  so-called  Minnesota  functionals  are,  in 
general,  not  as  accurate  as  D3,  but  they  can  perform  very  well  for  certain  systems  without 
the  need  for  any  additional  computational  effort.  This  is  particularly  true  for  compounds 
with  hydrogen-bonding  [56]— [59].  In  fact,  Grirnme  found  that  M06-2X,  a  hybrid 
Minnesota  functional,  performed  better  than  B3LYP-D3,  a  hybrid  functional  with  explicit 
D3  added,  when  benchmarking  hydrocarbon  reaction  energies,  which  include  chain,  ring 
and  cage  compounds  [60].  Explicit  inclusion  of  D3  with  M06-2X  only  produced  modest 
improvements  because  M06-2X  was  already  able  to  capture  most  of  the  dispersive  effects 
[60].  Grirnme  also  benchmarked  M06-2X  against  the  S66  and  S66x8  data  sets  which 
specifically  test  performance  for  noncovalent  interactions  [61].  He  found  excellent 
performance  using  M06-2X  for  H-bonded  complexes  with  very  little  improvement  with 
the  addition  of  D3  [61].  M06-2X  ran  into  difficulty  only  when  benchmarked  against  the 
S66x8  set  which  specifically  tests  long  range  dispersion  effects  by  separating  dimers 
beyond  their  equilibrium  geometry  [61].  Ultimately,  M06-2X  perfonns  well  for 
complexes  involving  short  and  medium  range  dispersion,  which  are  the  primary 
interactions  in  the  isolated  clusters  of  this  study.  It  is  import  to  point  out  that  the  success 
of  M06-2X,  or  any  other  functional,  will  be  highly  dependent  on  the  system  analyzed. 
This  is  also  true  of  the  techniques  that  are  designed  to  handle  dispersion.  While  explicit 
D3  will  more  often  than  not  provide  more  accurate  results  than  M06-2X,  this  is  not 
necessarily  the  case  for  every  system  [58],  In  fact,  recent  work  by  Lu  et  al.  has  shown 
that  M06-2X  has  proven  to  be  very  accurate  in  predicting  thermodynamic  properties  and 
key  structural  metrics  for  the  two  variants  with  reported  solid-state  structures  (the  Cp 
and  Cp  variants)  [62].  Because  of  this  and  the  computation  expediency  gained  by  not 
using  explicit  D3  corrections,  the  M06-2X  functional  was  chosen  as  the  primary 
workhorse  for  DFT  analysis. 

This  computational  approach  was  used  to  calculate  the  most  stable  isomers  of  the 
Cp  Pl  and  Cp  lPr  tetramers  in  the  absence  of  experimental  crystallographic  data.  Further 
details  on  determining  the  lowest  energy  geometric  configuration  are  discussed  in  II.A.4. 
Calculated  enthalpy  of  tetramerization  A Htet  values  for  these  isomers  are  within  3  kJ/mol 
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of  the  experimental  values  measured  for  the  two  new  structures.  Additionally,  M06-2X  is 
also  very  accurate  in  predicting  27A1  NMR  chemical  shifts  for  a  variety  of  low-valence 
aluminum  clusters.  Calculated  chemical  shifts  for  CpPl  and  Cp*lPr  tetramers  are  both 
within  9  parts  per  million  (ppm)  of  the  signals  measured  for  the  two  new  structures.  This 
strongly  suggests  that  the  two  new  structures  are  indeed  Cp  Pl  and  Cp*lPr  tetramers.  The 
M06-2X  functional  was  also  used  to  examine  the  thermodynamic  properties  for  the  other 
three  Cp-type  variants;  analysis  suggests  that  the  Cp  variant  has  a  significantly  lower 
tetramerization  enthalpy  compared  to  the  other  four  variants,  despite  previous  DFT  work 
showing  the  opposite  trend  [38]. 

A  detailed  analysis  of  bonding  was  conducted  in  the  context  of  the  Quantum 
Theory  of  Atoms  in  Molecules  (QTAIM)  theory.  This  analysis  shows  the  importance  of 
NCIs  between  the  ligands.  While  increased  ligand  bulk  can  result  in  increased  Al-Al 
separation  and  a  reduction  in  metal-metal  bonding  in  the  core,  this  is  often  offset  by  weak 
interactions  between  adjacent  ligand  groups.  The  tetramerization  enthalpies  arise  from  a 
balance  of  lone -pair  overlap  between  monovalent  aluminum  atoms  and  non-covalent 
ligand  interactions;  this  balance  is  challenging  for  DFT  functionals  to  capture  accurately. 
Finally,  using  thermodynamic  data  calculated  using  M06-2X,  monomer/tetramer 
equilibrium  constants  can  be  calculated  for  the  five  existing  homoleptic  Cp-type  variants 
as  well  as  five  theoretical  variants  with  varying  levels  of  steric  bulk. 
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II.  DFT  ANALYSIS  OF  MONOVALENT  ALUMINUM  CLUSTERS 


A.  TWO  NEW  C Y CLOPENT ADIEN YL  HOMOLEPTIC  CLUSTERS 

It  has  been  challenging  to  find  a  reliable  computational  method  for  monovalent 
aluminum  cluster  thermochemistry  using  ab  initio  methods,  due  in  large  part  to  a  lack  of 
experimental  data  and  the  difficult  synthesis  of  new  clusters  of  this  type.  Experimental 
structural  data  is  only  available  for  the  CpA  and  Cp  tetramers  and  the  larger  AfioCp  12, 
and  until  now  experimental  thermochemistry  data  was  only  available  for  the  Cp*  variant. 
Based  on  calculated  27A1  NMR  chemical  shifts  and  A Htet  values  for  AfiCp  '  '4  and 
AI4Cp  lPr4  it  is  assumed  that  experimental  values  observed  for  the  two  new  structures  are 
indeed  Cp*Pl  and  Cp  lPr  tetramers,  providing  two  additional  data  points  for  validating  the 
computational  accuracy.  Additionally,  these  two  compounds  are  important  because  of 
their  increased  ligand  steric  bulk  as  compared  to  the  Cp*variant.  Monomers  of  all  three 
variants  are  shown  in  Figure  9.  The  two  main  goals  of  this  thesis  are  to  demonstrate  the 
effectiveness  of  M06-2X  in  describing  the  properties  of  tetrameric  cyclopentadienyl 
aluminum  clusters  in  general  and  to  investigate  the  effect  of  varying  ligand  bulk  for  these 
clusters. 


Cp*  monomer  Cp*Pr  monomer  Cp*iPr  monomer 

The  aluminum  (blue)  is  positioned  on  a  line  perpendicular  from  the  middle  of  the  carbon  (grey)  ring. 


He  Hep.-  Hejp»- 

Figure  9.  Monomers  of  three  AIR  (R=Cp  ,  Cp  ,  Cp  )  variants. 
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1.  Synthesis  of  Al4Cp  Pl4  and  ALtCp*,Pl4 

The  synthesis  and  solution  equilibrium  of  the  C  p  Pr  and  Cp  lPr  variants  was  carried 
out  by  Dennis  Mayo  and  Rebecca  Wilson,  researchers  at  the  NSWC,  Indian  Head 
EODTD.  The  new  cyclopentadienyl  aluminum  derivatives  were  prepared  via  thermally- 
activated  reductive  elimination  of  tetramethyl(alkyl)cyclopentadiene  from  bis- 
(tetramethyl(R)cyclopentadienyl)  aluminum  hydride  intermediates  (Al(CsMe4R)2H,  R  = 
n-propyl  or  i-propyl)  in  toluene  at  80  °C  (see  Figure  10)  [63].  Higher  temperature 
reactions  were  attempted,  but  decomposition  products  were  observed  in  the  *H  NMR 
spectrum  when  carried  out  above  80  °C.  When  held  at  room  temperature,  solutions  of 
Al(CsMe4R)2H  in  toluene  or  benzene  are  also  in  equilibrium  as  observed  by  27A1  NMR 
spectroscopy:  a  solution  of  Al(C5Me4nPr)2H  had  an  Al(CsMe4nPr)2H:Al4Cp*Pl4  ratio  of 
65:35  after  two  weeks. 


Figure  10.  Synthesis  of  the  Cp  Pl  and  CplPr  variants  via  salt  metathesis  and 

reductive  elimination. 


Reductive  elimination  of  cyclopentadienes  from  Al(C5Me4R)2H  was  followed  via 
27A1  NMR  spectroscopy.  The  starting  Al(CsMe4R)2H  derivatives  have  diagnostic  27A1 
NMR  signals  around  -30  ppm;  the  Cp  Pl  and  Cp  lPr  oligomers,  assumed  to  be  tetramer 
forms,  have  27A1  signals  at  -86  and  -89  ppm,  respectively.  For  the  Cp  Pr  variant,  multiple 
cycles  of  heating  followed  by  removal  of  the  volatile  components  were  necessary.  For  the 
Cp  lPr  variant,  one  heat/evaporation  cycle  was  sufficient.  The  resultant  yellow  oils  and 
solutions  thereof  failed  to  yield  crystals  suitable  for  single  crystal  x-ray  analysis. 

The  resultant  cyclopentadienylaluminum  compounds  were  subjected  to  variable- 
temperature  van  ‘t  Hoff  analysis.  Solutions  of  the  Cp  Pl  and  Cp  lPr  variants  in  toluene -Js 
were  prepared  and  heated  to  80  °C,  the  27A1  and  *H  NMR  spectra  obtained,  and  the 
temperature  lowered  in  15  degree  increments.  Spectra  are  shown  in  Figure  11.  Upon 
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heating,  the  monomer  signal  grows  in  intensity  (-153  ppm  for  the  Cp1’1  variant,  -150  ppm 
for  the  Cp*lPr  variant),  showing  a  shift  in  the  equilibrium  towards  monomeric  species.  The 
thermodynamic  parameters  were  calculated  assuming  that  the  -86  and  -89  ppm  signals  for 
Cp  Pl  and  Cp  lPr  solutions  corresponded  to  tetrameric  forms;  attempts  to  fit  the  data 
assuming  monomer/dimer  equilibrium  resulted  in  considerably  worse  agreement. 
Additionally,  DFT  calculations  found  no  stable  dimer  or  trimer  configurations  of  AlCp  Pr 
or  AlCp*lPr.  A  dimer  configuration  of  AlCp*  that  was  artificially  frozen  in  place  was 
calculated  to  have  an  NMR  shift  of  -149  ppm,  almost  identical  to  that  of  the  monomeric 
fonn.  Thus,  while  larger  oligomers  cannot  be  ruled  out,  all  data  point  to  tetrameric  AI4R4 
structures.  For  the  Cp  Pr  variant,  the  enthalpy  and  entropy  changes  (A Htet  and  A Stet) 
derived  from  this  analysis  are  -160  ±3  kJ/mol  and  -398  ±10  J/mol-K,  respectively.  For 
the  Cp+lPr  variant,  A Htet  and  A Stet  are  -158  ±  8  kJ/mol  and  -All  ±  25  J/mol-K. 


Figure  11.  Variable-temperature  27A1  NMR  of  AfiCp  (a)  in  toluenc-rC  (100 
mg  in  0.8  mL;  302-347  K)  and  Al4Cp*lPr4  (b)  in  toluene-r/s  (44.5  mg  in 

0.5  mL;  303-353  K). 


2.  DFT  Approach  and  Thermochemistry  Results 

The  choice  of  DFT  functional  is  critical  regardless  of  the  type  of  system  under 
study.  Previous  theoretical  treatment  of  these  clusters  has  produced  a  range  of  different 
results.  DFT  calculations  by  Schnockel  and  co-workers  [38],  [41]  and  Williams  and 
Hooper  [31]  found  that  the  Cp  and  Cp  tetramers  are  both  much  more  stable  enthalpically 
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against  decomposition  into  monomers  than  the  Cp  tetramer.  This  seemed  to  be  supported 
by  the  experimental  fact  that  monomers  of  the  Cp  and  Cp  variants  were  not  found  via 
27A1  NMR  in  solution,  but  that  monomers  of  the  Cp  variant  were  found  [38],  [40].  While 
geometries  and  monomer  enthalpies  of  formation  were  accurately  predicted  in  these 
works,  tetramerization  enthalpies  were  in  poor  agreement  with  the  experimental  results 
available  for  the  Cp  variant.  Additionally,  new  experimental  thermochemistry  data  for 
the  Cp*Pl  and  Cp  'Pl  variants  showed  both  are  more  stable  enthalpically  than  the  Cp 
variant  even  though  monomers  of  the  Cp1’1  and  Cp  lPr  variants  are  found  in  solution. 
With  this  new  experimental  thermochemistry  data  available,  an  examination  of  the  Cp*, 
Cp  Pl  and  Cp  lPr  variants  was  conducted  with  a  range  of  functionals  and  basis  sets  to 
determine  a  computational  approach  that  is  accurate  for  both  structural  and 
thermochemical  values.  Since  crystal  structures  do  not  exist  for  the  Cp*Pr  and  Cp  lPr 
tetramers,  multiple  configurations  of  each  were  examined  to  determine  the  global 
minimum  cluster  structure  and  nearby  isomers.  That  process  is  discussed  later  in  this 
work. 

Ten  functionals  [55],  [64]-[74]  were  examined,  all  with  a  split  valence  double- 
zeta  basis  sets  with  polarization.  This  type  of  basis  set  was  used  for  computational 
expediency.  Included  in  this  list  are  the  two  methods  used  in  previous  works  on  these 
clusters.  The  first  is  the  method  of  Schnockel  which  is  the  GGA  functional  BP86  paired 
with  the  split  valence  double-zeta  basis  set  with  polarization,  SVP  [38].  The  second  is  the 
method  of  Williams  and  Hooper  which  is  the  hybrid  functional  B3LYP  paired  with  the 
split  valence  double-zeta  basis  set  with  polarization,  6-31G(d,p)  [31].  All  functionals 
produce  accurate  results  for  key  geometrical  parameters  such  as  the  Al-Al  distance  and 
the  Al-Cp  ring-center  distance.  Structural  parameters  for  all  functionals  are  given  in 
Table  1  and  in  general  these  predict  key  distances  to  better  than  0.05  A  compared  to 
solid-state  x-ray  data  for  the  Cp*  tetramer.  To  get  a  clearer  picture  of  overall  functional 
performance,  Table  2  shows  the  differences  in  calculated  geometrical  parameters  from 
experimental  ones.  While  all  functionals  perform  well,  M06-2X  is  the  most  accurate  for 
three  of  the  five  metrics  and  is  the  second  and  third  most  accurate  for  the  remaining  two. 
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It  is  also  important  to  note  that  M06-2X  performs  better  than  B3LYPD3,  a  hybrid 
functional  with  explicit  dispersion  corrections  added,  in  all  five  metrics. 


Table  1.  Calculated  average  lengths  for  AfiCpV 


Method 

Al-Al 

Al-CpCenter 

Al-Cring 

C  ring- C  ring 

Cring“CMe 

EXPERIMENTAL 

2.7670 

2.0110 

2.3450 

1.4170 

1.5000 

M062X/6-31G(d) 

2.7672 

1.9925 

2.3313 

1.4231 

1.5003 

BP86/SVP 

2.8184 

2.0570 

2.3935 

1.4386 

1.5041 

B3LYP/6-3  lG(d,p) 

2.8541 

2.0693 

2.3994 

1.4278 

1.5029 

B3LYPD3/6-3  lG(d,p) 

2.7403 

2.0295 

2.3646 

1.4268 

1.5018 

TPSSh/6-31G(d) 

2.7802 

2.0495 

2.3830 

1.4293 

1.5037 

PBEPBE/6-3  lG(d) 

2.7989 

2.0505 

2.3865 

1.4355 

1.5030 

PBE 1 PBE/6-3 1  G(d) 

2.7849 

2.0273 

2.3615 

1.4239 

1.4957 

M06L/6-31G(d) 

2.6765 

1.9801 

2.3214 

1.4248 

1.4920 

mPWlPW9 1/6-3 1  G(d) 

2.7989 

2.0367 

2.3694 

1.4235 

1.4962 

APFD/6-31G(d) 

2.7234 

1.9988 

2.3370 

1.4237 

1.4968 

*  All  numbers  are  in  angstroms. 


Table  2.  Calculated  length  differences  from  experiment  for  AfiCpV 


Method 

Al-Al 

Al-CpCenter 

Al-Cring 

C  ring- C  ring 

Cring-C]y[e 

M062X/6-31G(d) 

0.0002 

0.0185 

0.0137 

0.0061 

0.0003 

BP86/SVP 

0.0514 

0.0460 

0.0485 

0.0216 

0.0041 

B3LYP/6-3  lG(d,p) 

0.0871 

0.0583 

0.0544 

0.0108 

0.0029 

B3LYPD3/6-3  lG(d,p) 

0.0267 

0.0185 

0.0196 

0.0098 

0.0018 

TPSSh/6-31G(d) 

0.0132 

0.0385 

0.0380 

0.0123 

0.0037 

PBEPBE/6-3  lG(d) 

0.0319 

0.0395 

0.0415 

0.0185 

0.0030 

PBE  1  PBE/6-3  lG(d) 

0.0179 

0.0163 

0.0165 

0.0069 

0.0043 

M06L/6-31G(d) 

0.0905 

0.0309 

0.0236 

0.0078 

0.0080 

mP  W 1 P  W  91/6-3  lG(d) 

0.0319 

0.0257 

0.0244 

0.0065 

0.0038 

APFD/6-31G(d) 

0.0436 

0.0122 

0.0080 

0.0067 

0.0032 

*  All  numbers  are  in  angstroms. 


Despite  accurate  geometrical  results  for  all  functionals,  thermodynamic  results  are 
much  more  sensitive  to  functional.  Figure  12  shows  the  strong  variation  of  the 
tetramerization  enthalpy  (A Htet)  with  functional  type.  All  three  experimental  clusters 
were  optimized  with  each  functional  until  no  imaginary  frequencies  were  found.  The 
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enthalpy  and  entropy  for  the  tetramerization  reaction  were  calculated  from  vibrational 
frequency  analysis. 


2<K) 


BPH6  B3LYP  TPSSh  PBEPBE  PBEO  M06L  mPWlPW9l  A  PHD  M062X 

SVP  6-3IG(d,p)  6-3IG<d>  6-3IG(d)  6-3IG(d)  6-31G(d)  6-3IG(dt  A-3IG(d)  6-3IG(d) 

Method 


Figure  12.  Difference  of  DFT  calculated  AHtet  values  from  experimental  values 

for  various  functionals. 


Differences  in  A Htet  between  experiment  and  theory  are  as  high  as  130  kJ/mol,  a 
value  sufficient  to  make  equilibrium  rate  calculations  impossible.  Additionally,  for  all  but 
three  functionals,  the  ordering  of  A Htet  for  the  three  variants  differ  from  experiment.  The 
Minnesota  LDA  functional,  M06-L,  and  the  hybrid  functional,  APFD,  produce  the 
correct  relative  ordering  of  tetramers,  but  both  dramatically  over-predict  A Htet  for  all 
three  structures.  The  third,  M06-2X,  has  results  much  closer  to  experiment.  While  PBE, 
PBEO  and  M06-2X  all  perform  significantly  better  than  the  other  methods  tested,  only 
M06-2X  predicts  the  correct  order  of  A Htet  for  all  three  variants.  Capturing  this  trend  is 
what  is  important,  as  the  double-zeta  basis  set  is  not  expected  to  be  of  high  accuracy.  The 
results  for  A Stet  are  more  mixed.  Two  methods,  BP86  and  B3LYP  get  the  order  correct, 
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but  the  calculated  result  for  the  Cp  variant  is  in  error  by  over  100  J/mol-K.  In  general,  it 
appears  all  tested  methods  have  much  higher  errors  for  predicting  A Stet.  Numerical 
results  for  A Htet  and  A Stet  are  shown  in  Table  3.  In  detailed  analysis  to  follow,  enthalpy 
changes  are  relied  on  to  validate  predicted  structures. 


Table  3 .  Heats  and  entropies  of  tetramization  for  the  Cp  ,  Cp  Pl  and  Cp*lPl 
variants  using  various  DFT  functionals. 


Method 

AfiCp  4 

Al4Cp*Pr4 

Al4Cp*iPr4 

A  Htet 

A  Htet 

AStet 

A  Htet 

A  Stet 

EXPERIMENTAL 

-150 

-0.300 

-160 

-0.398 

-158 

-0.477 

BP86/SVP 

-100 

-0.456 

-98 

-0.482 

-82 

-0.563 

B3LYP/6-3  lG(d,p) 

-17 

-0.436 

-27 

-0.453 

-10 

-0.585 

TPSSh/bGIGtd)1 

-126 

-0.523 

-122 

-0.496 

-103 

-0.577 

PBEPBE/6-3  lG(d) 

-163 

-0.524 

-164 

-0.503 

-150 

-0.537 

PBE 1 PBE/6-3 1  G(d) 

-150 

-0.497 

-150 

-0.497 

-137 

-0.560 

M06-L/6-31G(d) 

-285 

-0.564 

-301 

-0.533 

-293 

-0.665 

mP  W 1 P  W  9 1/6-3 1  G(d) 

-117 

-0.526 

-116 

-0.514 

-100 

-0.544 

APFD/6-31G(d) 

-241 

-0.532 

-256 

-0.613 

-255 

-0.633 

M06-2X/6-3  lG(d) 

-174 

-0.649 

-188 

-0.679 

-186 

-0.657 

*  A Htet  numbers  are  in  kJ/mol;  AStet  numbers  are  in  kJ/mol-K 


To  determine  the  sensitivity  of  basis  set  choice  on  thermodynamic  calculations, 
the  M06-2X  functional  was  paired  with  eight  different  basis  sets  from  Pople,  Dunning 
and  Ahlrichs  [75]— [80].  The  resulting  thermodynamic  calculations  for  A Htet  and  A Stet,  in 
addition  to  experimental  results  for  each  of  the  three  variants,  are  shown  in  Table  4. 
These  results  do  not  include  two  energy  corrections  required  for  these  clusters.  The  first 
is  basis  set  super  position  error  (BSSE).  BSSE  corrects  for  artifacts  that  arise  when 
molecular  energy  is  calculated  for  an  isolated  monomer  and  a  monomer  bound  in  a 
cluster.  The  second  is  energy  corrections  due  to  the  presence  of  the  toluene  solvent.  Both 


1  Despite  multiple  techniques,  the  last  imaginary  vibrational  frequency  for  the  Cp*  variant’s  tetramer 
using  TPSSh  was  not  eliminated.  Based  on  experience  reducing  negative  vibrational  frequencies  from  all 
other  methods  it  was  determined  that  the  removal  of  a  negative  frequency  corresponds  to  roughly  a  2.5 
kJ/mol  reduction  in  A Htet  and  about  a  0.150  kJ/mol-K  reduction  in  A Stet.  These  corrections  are  not 
reflected  in  Table  3,  but  would  not  change  our  overall  analysis  of  functional  performance. 
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corrections  were  omitted  for  computational  expediency.  While  there  is  not  a  systematic 
approach  to  accurate  A Htet  values  with  basis  set  size,  the  larger  cc-PVTZ  was  chosen  as 
the  most  accurate  available  method.  Calculations  of  tetramers  at  the  double  zeta  level 
typically  required  approximately  4000  CPU  hours  per  system  on  a  Department  of 
Defense  (DOD)  High  Perfonnance  Computing  (HPC)  cluster;  increasing  the  basis  size  to 
cc-PVTZ  required  up  to  25,000  CPU  hours  on  the  same  machine.  Calculations  on  larger 
clusters,  such  as  AfioCp  12.  were  not  viable  at  the  triple  zeta  level  with  the  available 
resources  on  DOD  HPC  systems.  For  this  reason,  and  due  to  the  ability  of  M06-2X/cc- 
pVTZ  to  predict  A Htet  to  within  3  kJ/mol  of  experimental  values  for  all  variants  when 
BSSE  and  solvent  corrections  are  added,  calculations  beyond  the  triple  zeta  level  of 
theory  were  not  perfonned. 


Table  4.  Heats  and  entropies  of  tetramization  for  the  Cp  ,  Cp  Pl  and  Cp*lPl 
variants  using  M06-2X  and  various  basis  sets. 


Method 

AI4Cp*4 

AI4Cp*Pr4 

AI4Cp*iPr4 

A  Htet 

h^tet 

A  Htet 

h^tet 

A  Htet 

A  Stet 

EXPERIMENTAL 

-150 

-0.300 

-160 

-0.398 

-158 

-0.477 

M062X/6-31G(d) 

-174 

-0.649 

-188 

-0.679 

-186 

-0.657 

M062X/6-31+G(d) 

-168 

-0.624 

-182 

-0.678 

-180 

-0.679 

M062X/6-31G(d,p) 

-181 

-0.606 

-193 

-0.676 

-189 

-0.648 

M062X/6-311G(d,p) 

-172 

-0.650 

-188 

-0.626 

-187 

-0.653 

M062X/6-311+G(d,p) 

-172 

-0.658 

-185 

-0.635 

-183 

-0.650 

M062X/Def2-TVZP 

-173 

-0.633 

-183 

-0.690 

-179 

-0.671 

M062X/cc-pVDZ 

-168 

-0.610 

-181 

-0.666 

-180 

-0.658 

M062X/cc-pVTZ 

-167 

-0.639 

-178 

-0.690 

-175 

-0.663 

*  A Htet  numbers  are  in  kJ/mol;  A Stet  numbers  are  in  kJ/mol-K 


Table  5  shows  calculated  A Htet  values  when  D3  is  included  with  M06-2X  and 
B3LYP.  The  addition  of  Grimme’s  D3  dispersion  correction  decreases  the  predicted 
A Htet  values  calculated  with  M06-2X  by  ~40  kJ/mol.  This  is  not  surprising  as  M06-2X 
has  already  been  specifically  parameterized  to  capture  NCIs.  Including  D3  with  B3LYP 
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dramatically  improved  results,  but  the  functional  still  incorrectly  orders  the  Cp  Pl  and 

_ ,  *iPr 

Cp  variant. 


Table  5.  Heats  of  tetramization  for  the  Cp*,  Cp*Pr  and  Cp*lPr  variants  with 

added  dispersion. 


Method 

AfiCp  4 

Al4Cp*Pr4 

Al4Cp*lPr4 

EXPERIMENTAL 

-150 

-160 

-158 

B3LYPD3/6-3  lG(d,p) 

-167 

-180 

-184 

M062XD3/6-3  lG(d) 

-202 

-223 

-220 

*  A Htet  numbers  are  in  kJ/mol. 


The  M06-2X  functional  is  clearly  superior  to  all  other  methods  tested  and  does 
not  require  explicit  treatment  for  van  der  Waals  interactions.  This  choice  is  consistent 
with  recent  work  done  on  the  Cp*  and  Cp  variants  as  well  as  the  related  heterocubanes 
reported  by  Stelzer  and  co-workers  [45],  [62],  Table  6  shows  the  final  calculated 
thermodynamic  values  for  all  five  clusters  using  M06-2X/cc-pVTZ  including  BSSE  and 
solvent  corrections. 


Table  6.  Calculated  A Htet  including  basis  set  superposition  error  correction 

for  all  variants. 


AfiCp  4 

Al4Cp4 

AfiCp  4 

Al4Cp*Pr4 

Al4Cp*lPr4 

M062X/cc-pVTZ 

-148.4 

-106.7 

-153.1 

-158.7 

-155.8 

Experimental 

-150 

- 

- 

-160 

-158 

*  All  numbers  are  in  kJ/mol. 


3.  NMR  Results 

In  addition  to  accurate  enthalpies  of  tetramerization,  predicted  27A1  NMR  shifts 
using  the  M06-2X  functional  are  quite  accurate  for  known  monovalent  aluminum 
clusters.  There  are  eight  cyclopentadienyl  aluminum  structures  discussed  in  this  work 
with  experimental  27A1  NMR  data:  the  Cp  and  Cp  tetramers  and  the  monomers  and 
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tetramers  of  the  Cp*,  Cp  Pl  and  Cp*lPl  variants.  In  addition  to  these  eight,  experimental 
NMR  data  for  nine  additional  related  aluminum  cyclopentadiene  clusters  have  been 
reported  by  Sitzmann  and  coworkers  [35].  Gauge-independent  atomic  orbital  (GIAO) 
calculations  with  M06-2X  using  a  6-31G(d)  basis  set  were  performed  and  the  results 
show  that  a  static  shift  of  +4.5  ppm  brings  calculated  results  into  close  agreement  with 
experiment.  Test  calculations  with  the  larger  cc-pVTZ  basis  set  using  the  M06-2X 
functional  showed  virtually  the  same  results,  just  with  a  slightly  larger  static  shift.  A  plot 
of  experimental  values  against  calculated  values,  and  associated  regression  lit,  are  shown 
in  Figure  13. 


Figure  13.  Calculated  vs.  experimental  chemical  shifts  using  M06-2X/6-31G(d). 

The  use  of  other  functionals  showed  strong  variations  from  experiment.  Figure 
14,  for  example,  shows  the  poor  predictive  capability  of  a  functional  such  as  BP86  which 
was  used  in  previous  works.  The  plot  shows  the  same  17  structures  from  Figure  13,  but 
the  calculated  ppm  values  were  determined  with  BP86/SVP.  A  linear  regression  results  in 


30 


a  y-intercept  of  -35,  but  also  in  a  slope  of  0.7.  This  slope  indicates  a  systematic  error 
present  when  using  the  BP86/SVP  method.  Specifically,  at  smaller  ppm  magnitudes, 
calculated  results  tend  to  be  more  negative  than  experimental  results.  At  about  -115  ppm, 
calculated  and  experimental  tend  to  match.  At  magnitudes  beyond  -115  ppm,  calculated 
results  tend  to  be  less  negative  than  experimental  results.  Additionally,  the  linear  fit  with 
an  R2  of  0.9286  is  not  as  good  as  the  one  for  M06-2X  with  an  R2  of  0.978 1 . 


Experiment  5  (ppm) 

Figure  14.  Calculated  chemical  shifts  using  BP86/SVP. 

To  illustrate  this  point,  Table  7  shows  the  calculated  and  experimental  NMR 
chemical  shifts  for  all  17  structures.  It  is  ordered  from  the  most  negative  experimental 
value  to  least  negative  and  includes  the  difference  between  experiment  and  calculation 
for  each  structure.  For  BP86/SVP,  the  average  difference  from  experiment  for  the  seven 
structures  with  chemical  shifts  less  than  -115  ppm  is  17.5.  The  average  difference  from 
experiment  for  the  ten  structures  with  chemical  shifts  greater  than  - 1 1 5  ppm,  however,  is 
-10.1.  For  M06-2X/6-31G(d),  the  average  difference  from  experiment  for  the  seven 
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structures  with  chemical  shift  less  than  -115  ppm  is  8.5.  For  the  ten  structures  with 
chemical  shifts  more  than  -115  ppm  the  average  is  4.3.  This  is  consistent  with  assertion 
that  M06-2X/6-31G(d)  provides  NMR  calculations  that  are  generally  about  5  ppm  off 
from  experiment. 


Table  7.  Difference  of  calculated  7A1  NMR  chemical  shifts  from  experiment. 


Structure 

Calculated  Results  with  de 

Ita 

Experiment 
al  Results 

BP86 

delta 

M062X 

delta 

Al(C5H3(SiMe3)2) 

-151.0 

17.0 

-157.4 

10.6 

-168 

Al(C5H2(SiMe3)3) 

-147.3 

17.7 

-152.8 

12.2 

-165 

Al(C5/Pr4H) 

-135.6 

23.4 

-144.8 

14.2 

-159 

AlC5(CH2Ph)5 

-137.9 

17.1 

-154.9 

0.1 

-155 

* 

Cp  monomer 

-135.3 

14.7 

-144.2 

5.8 

-150 

p  *iPr 

Cp  monomer 

-136.3 

13.7 

-141.6 

8.4 

-150 

Cp  Pl  monomer 

-134.4 

18.6 

-144.5 

8.5 

-153 

Cp  tetramer 

-117.3 

-6.3 

-110.8 

0.2 

-111 

(AlCp*)3AlCps 

-123.7 

-15.7 

-104.5 

3.5 

-108 

Al4(C5H4fflu)4 

-118.4 

-11.4 

-107.7 

-0.7 

-107 

Al4(C5H3(SiMe3)2)4 

-119.7 

-13.7 

-109.3 

-3.3 

-106 

Cp  tetramer 

-93.0 

-4.0 

-90.8 

-1.8 

-89 

Cp  tetramer 

-99.4 

-13.4 

-74.8 

11.2 

-86 

Cp  tetramer 

-93.5 

-11.5 

-74.2 

7.8 

-82 

Cp  tetramer 

-92 

-11 

-71.3 

9.7 

-81 

(AlCp*)3AlCpCT 

-82.6 

-5.6 

-68.2 

8.8 

-77 

(AlCp*)3(AlN(SiMe3)2) 

-70.3 

-8.3 

-55.4 

6.6 

-62 

s  A1  attached  to  Cp  ligand 


s$  A1  attached  to  Cp*  ligand 


4.  Structural  Analysis  of  Al4Cp  Pr4  and  Al4Cp  lPr4 

Solutions  the  Cp  Pl  and  Cp  lPr  variants  failed  to  produce  x-ray  quality  crystals,  and 
thus  direct  structural  infonnation  is  not  available  for  these  two  newly  synthesized 
systems.  Therefore,  multiple  calculations  at  the  M06-2X/cc-pVTZ  DFT  level  of  theory 
were  completed  and  resulting  tetramerization  energies  were  compared  with  variable 
temperature  NMR  experiments.  For  computational  expediency,  solvent  corrections  were 
not  calculated. 
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The  additional  ligand  steric  bulk  in  the  Cp  Pl  and  Cp  lPr  variants  results  in  a 
number  of  potential  isomeric  configurations  of  the  tetramers.  Unlike  the  Cp  monomers, 
the  monomers  of  the  Cp  Pl  and  Cp  'Pl  variants  are  not  symmetric.  DFT  is  therefore  used  to 
determine  the  ground  state  structure,  beginning  with  the  Cp1’1  variant.  Two  stable 
configurations  of  the  monomeric  AlCp  Pr  were  found  and  are  displayed  in  Figure  15.  In 
the  first,  shown  in  Figure  15(a)  and  referred  to  here  as  type  A,  the  propyl  group  is 
oriented  away  from  the  monovalent  aluminum.  The  type  B  configuration  shown  in  Figure 
15(b)  has  the  propyl  directed  inwards  in  a  configuration  that  would  be  expected  to  lead  to 
increased  steric  hindrance  in  the  tetramer. 


The  type  A  monomer  is  only  1.9  kJ/mol  lower  in  DFT  energy  than  the  type  B 
monomer.  Both  types  of  monomers  were  found  to  combine  into  stable  tetramers,  though 
the  tetramer  formed  by  type  B  monomers  is  less  enthalpically  stable  than  its  type  A 
counterpart.  The  most  enthalpically  stable  tetramers  derived  from  each  variant  are  shown 
in  Figure  16.  The  type  B  tetramer  has  a  BSSE  corrected  A Htet  of  -135.3  kJ/mol  as 
compared  to  the  -165.4  kJ/mol  for  the  type  A  tetramer  without  the  solvent  correction. 
This  appears  to  be  primarily  due  to  a  larger  Al-Al  separation  in  the  type  B  tetramer,  as 
might  be  expected.  The  average  Al-Al  distance  in  the  type  B  variant  is  2.867  A,  10 
picometers  larger  than  the  2.767  A  distance  for  type  A  tetramers,  and  is  consistent  with 
weaker  bonding.  Also  considered  were  tetrameric  variants  formed  with  a  mix  of  type  A 
and  B  monomers.  While  none  of  these  were  found  to  have  the  highest  enthalpic  stability, 
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their  predicted  A Htet  values  are  close  to  the  most  enthalpically  stable  structure  tested  and 
the  experimental  value.  The  most  enthalpically  stable  mixed  configuration,  a  tetramer 
composed  of  one  B  type  monomer  and  three  A  type  monomers,  has  a  A Htet  of  -161.9 
kJ/mol,  just  over  3  kJ/mol  lower  than  a  tetramer  composed  of  all  A  type  monomers.  It 
also  has  an  average  Al-Al  distance  of  2.794  A  which  is  only  slightly  larger  than  the 
average  Al-Al  distance  in  the  pure  type  A  tetramer. 


Hydrogen  atoms  removed  for  clarity. 

Figure  16.  Type  A  (a)  and  type  B  (b)  tetramers  for  the  Cp  Pr  variant. 

Beyond  the  difference  in  monomers  of  the  Cp  Pr  variant,  there  are  also  many 
possible  configurations  of  the  Cp  Pr  tetramer.  To  understand  these  configurations,  it  is 
helpful  to  imagine  the  tetramer  as  a  tetrahedron,  with  each  face  of  the  tetrahedron 
represented  by  a  Cp  Pl  ligand.  There  are  two  distinct  directions  the  propyl  arm  of  the 
ligand  can  point.  It  can  point  towards  the  edge  of  an  adjacent  ligand,  in  which  case  the 
propyl  arm  settles  between  two  adjacent  methyl  groups  from  the  neighboring  ligand,  or  it 
can  point  towards  the  gap  between  two  adjacent  ligands,  which  is  essentially  one  of  the 
vertexes  of  the  tetrahedron.  Each  type  of  configuration,  edge  or  vertex,  has  separate 
configurations  within  it.  The  configurations  are  visualized  in  Figure  17  by  a  flattened 
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tetrahedron  representing  the  four  faces  of  the  central  AI4  tetrahedron  and  their  respective 
Cp  Pr  configurations.  DFT  energies  are  listed  for  each  configuration  relative  to  the  global 
minimum. 


(e)  +  7.2  kJ/mol 


(b)  +4.3  kJ/mol 


(d)  +6.1  kJ/mol 


(g)  +9.3  kJ/mol 


(h)  +11.8  kJ/mol 


Me  groups  omitted  for  clarity. 

Figure  17.  Eight  different  isomers  of  the  Cp  Pr  tetramer,  ordered  by  DFT 

energies. 


The  A Htet  of  these  configurations  are  all  within  13  kJ/mol  of  each  other,  which  is 
comparable  to  the  error  in  DFT  thermochemical  results  (3  kJ/mol).  The  A Htet  of  the  most 
enthalpically  stable  configuration  shown  in  Figure  17(a),  however,  is  nearly  5  kJ/mol 
lower  than  the  next  closest  isomer,  Figure  17(b). 

Unlike  the  Cp*Pl  monomer,  Cp*lPl  monomers  are  symmetric  about  the  plane  of  the 
ligand  and  do  not  have  multiple  variants.  Calculations  of  isomers  showed  that  Cp  lPr 
configurations  in  which  the  iPr  group  points  towards  the  vertex  of  the  tetrahedron  are 
strongly  preferred.  Because  part  of  the  isopropyl  group  protrudes  towards  the  aluminum 
core,  arranging  four  ligands  in  one  of  the  various  edge  configurations  forces  a  slight 
separation  between  the  four  core  aluminum  atoms,  which  results  in  a  significant  increase 
in  the  calculated  A Htet  that  is  inconsistent  with  experiment.  A  variant  similar  to  Figure 
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17(a),  the  most  enthalpically  stable  configuration  for  Cp  Pr,  resulted  in  a  predicted  A Htet 
of -135.6  kJ/mol,  significantly  different  than  the  experimental  value  of -158  kJ/mol.  The 
Al-Al  bond  distance  for  this  configuration  was  found  to  be  2.863  A,  almost  0.1  A  further 
apart  than  the  Al-Al  bond  distance  of  the  Cp  Pl  variant  in  this  configuration.  The  A Htet 
for  vertex  configurations  were  much  closer  to  experimental  results.  In  addition  to 
analyzing  configurations  (c),  (f),  (g)  and  (h)  from  Figure  17  for  the  Cp  lPr  variant,  one 
additional  vertex  configuration  was  analyzed.  This  configuration  is  very  similar  to  (f) 
except  that  instead  of  two  sets  of  two  isopropyl  arms  pointing  to  the  same  vertex  there  is 
only  one  set  of  two  isopropyl  arms  pointing  towards  the  same  vertex.  The  other  two 
isopropyl  arms  point  to  separate  vertices.  This  configuration,  (i),  is  shown  in  Figure  18. 
The  calculated  values  of  A Htet  for  the  five  configurations  were  all  within  8  kJ/mol  of 
each  other,  with  configuration  (i)  the  most  enthalpically  stable.  Of  note,  the  Al-Al  bond 
distances  for  this  configuration  was  2.824  A.  While  this  distance  is  still  greater  than  the 
Al-Al  bond  distance  of  the  Cp  Pl  variant,  it  is  about  40  percent  closer  than  the  isopropyl 
edge  configuration  analyzed. 


(') 

Me  groups  omitted  for  clarity. 

Figure  18.  The  lowest  energy  configuration  for  the  Cp*lPl  variant. 
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B.  ANALYSIS  OF  THE  CP*ET  VARIANT  AND  OVERALL  LIGAND  TRENDS 

Using  this  same  approach,  another  cluster  in  this  monovalent  aluminum  series, 
Al4Cp  Et4  (Cp  Et=C5Me4Et),  was  examined.  While  this  system  has  yet  to  be  synthesized,  it 
allows  for  the  analysis  of  trends  in  the  tetramer  properties  as  a  function  of  Cp  bulk.  The 
most  enthalpically  stable  configuration  for  this  variant  is  shown  in  Figure  19.  This  is  the 
same  configuration  as  the  lowest  energy  Cp*Pl  variant  shown  in  Figure  17(a).  The  A Htet 
for  this  cluster  is  -157.9  kJ/mol,  which  makes  it  slightly  more  strongly  bound  than  the 
Cp*  variant,  but  less  strongly  bound  than  the  CpA  variant. 


Me  groups  omitted  for  clarity. 

Figure  19.  The  lowest  energy  configuration  for  the  theoretical  Cp11  variant. 

After  structural  analysis  of  the  Cp  Pr,  Cp*lPl  and  Cp  Et  variants  was  complete,  all 
six  clusters  discussed  up  to  this  point  were  optimized  with  M06-2X/cc-pVTZ.  BSSE 
corrections  were  added  as  well  as  corrections  based  on  the  presence  of  the  solvent 
toluene.  Solvent  corrections  increased  A Htet  values  by  about  6  kJ/mol  for  all  clusters,  but 
brought  calculated  values  very  close  to  experimental  results.  Figure  20  shows  the  most 
enthalpically  stable  structures  for  all  six  Cp  variants,  listed  in  order  of  increasing  steric 
bulk  and  annotated  with  calculated  A Htet  values.  Experimental  values  are  also  included 
for  the  Cp*,  Cp*Pl  and  Cp*lPr  variants.  In  general,  adding  steric  bulk  within  this  series 
results  in  tetramers  with  stronger  binding,  but  AfiCp  4  is  an  exception  to  this  trend. 
Crystal  structure  data  for  the  Cp  variant  show  that  it  is  configured  as  shown  in  Figure 
17(b),  though  it  should  be  noted  that  DFT  calculations  show  that  the  Cp  variant 
configured  as  in  Figure  17(a)  has  a  A Htet  only  0.5  kJ/mol  higher  than  the  crystal 


37 


structure.  This  suggests  that  the  Cp  variant  may  be  somewhat  insensitive  to 
configuration.  Despite  the  fact  that  the  Cp*  and  Cp*Et  variants  carry  more  steric  bulk,  the 
Cp  variant  is  more  enthalpically  stable.  Stability  gained  by  inter-ligand  interactions  can 
be  offset  by  a  loss  of  Al-Al  bond  strength  due  to  increased  monomer  separation  within  a 
cluster.  To  better  understand  this  balance,  the  electronic  structure  is  considered  in  more 
detail. 


AI4CP4 

calc:  -106.7  kJ/mol 
exp:  none 


calc:  -153.1  kJ/mol 
exp:  none 


AI4Cp  4 

calc:  -148.4  kJ/mol 
exp:  -150  kJ/mol 


AI4Cp‘Et4 

calc:  -151.5  kJ/mol 
exp:  none 


calc:  -155.8  kJ/mol 
exp:  -158  kJ/mol 


calc:- 158.7  kJ/mol 
exp:  -160  kJ/mol 


Calculated  and  experimental  A Htet  are  included.  Hydrogens  removed  for  clarity. 


Figure  20.  Initial  six  variants  investigated  in  order  of  increasing  steric  bulk. 
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C.  QUANTUM  THEORY  OF  ATOMS  IN  MOLECULES  ANALYSIS 

While  Al-Al  bond  length  within  the  core  of  these  clusters  can  inform  overall  bond 
strength,  this  length  does  not  appear  to  be  good  metric  for  gauging  overall  cluster 
stability.  The  lowest  energy  configuration  of  the  Cp  lPr  variant,  for  example,  has  larger 
Al-Al  separation  than  the  Cp,  Cp  and  Cp*  variants  but  is  enthalpically  preferred.  Work 
done  by  Stelzer  and  co-workers  demonstrated  the  existence  of  attractive  NCIs  between 
adjacent  ligands  in  the  Cp  variant  [45].  Additionally,  Lu  and  co-workers  demonstrated 
the  existence  of  attractive  NCIs  across  ligands  in  the  Cp*  variant  [62],  To  understand  the 
role  of  inter-ligand  bonding  in  all  variants,  QTAIM  analysis  is  performed. 

QTAIM  analysis  involves  examining  key  topological  parameters  of  the  electronic 
density  of  a  given  molecule.  In  particular,  points  where  the  density  gradient  vanishes  can 
indicate  potential  bonds  between  atoms.  Figure  21  shows  a  plot  of  the  norm  of  the 
gradient  for  the  Cp"  anion.  All  values  are  positive  and  are  color  coded  such  that  red  is  a 
maximum  and  blue  is  zero.  Within  the  Cp"  molecule  there  are  1 1  critical  points  where  Vp 
vanishes.  Ten  of  these  critical  points  clearly  correspond  with  the  ten  covalent  bonds 
present  in  the  molecule.  These  critical  points  are  three-dimensional  saddle  points  with 
values  of  p  decreasing  in  two  orthogonal  directions  and  increasing  in  the  third.  These 
critical  points  are  called  bond  critical  points  (BCPs)  because  they  correspond  to 
maximum  pathways  of  electron  density  between  two  atoms.  The  1 1th  critical  point,  at  the 
center  of  the  molecule,  is  a  three-dimensional  saddle  point  with  values  of  p  increasing  in 
two  orthogonal  directions  and  decreasing  in  the  third.  This  critical  point  is  called  a  ring 
critical  point  and  is  purely  consequence  of  the  ring  structure  of  the  molecule.  Using  these 
and  other  metrics  of  p,  QTAIM  provides  a  means  of  locating  bonds  and  a  quantitative 
way  to  classify  them.  Of  particular  importance  is  the  ability  of  QTAIM  to  predict  NCIs 
that  are  difficult  if  not  impossible  to  predict  from  molecular  structure.  These  types  of 
interactions  are  key  contributors  to  the  overall  stability  of  the  clusters  studied  in  this 
work. 
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Figure  2 1 .  Plot  of  the  norm  of  Vp  for  the  Cp"  anion. 

BCPs  indicating  NCIs  are  found  across  ligands  in  the  Cp  variant  and  values  for 
electron  density  at  these  and  other  key  BCPs  ( pb )  are  within  one  percent  of  previous 
work  [62].  Similar  BCPs  are  found  across  ligands  of  the  two  newly  synthesized  Cp  Pl  and 
Cp  lPr  tetramers.  Figures  22  and  23  show  BCPs  between  adjacent  ligands  for  the  lowest 
energy  tetramers  for  both  variants.  BCPs  between  hydrogen  atoms  as  well  as  hydrogen 
and  carbon  atoms  across  ligands  are  observed.  No  inter-ligand  bonding  was  found  for  the 
Cp  variant.  BCPs  between  hydrogen  atoms  within  a  single  ligand  are  also  observed,  but 
these  do  not  contribute  to  overall  tetramer  stability. 


Figure  22.  Inter-ligand  H-H  BCPs  in  the  Cp  Pl  variant  (shown  in  red). 
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Figure  23.  Inter-ligand  H-C  and  H-H  BCPs  in  the  Cp  'Pl  variant  (shown  in  red). 

All  inter-ligand  BCPs  were  characterized  by  low  values  of  pb,  positive  values  of 
V2pb  and  a  positive  value  for  total  energy  density  at  the  BCP,  Hh.  These  are  all  indicators 
of  closed  shell  bonding,  or  more  generally,  NCIs  [81],  [82].  The  QTAIM  parameters  for 
Al-Al  BCPs  indicate  a  shared  or  covalent  type  bond.  This  is  due  to  V2pb  <  0  and,  Hb  < 
0,  both  indicators  of  a  covalent  bond  [83],  [84].  The  fact  that  pb  at  these  points  is  slightly 
lower  than  the  more  typical  10"1  e-A"3  only  indicates  that  the  shared  bond  is  relatively 
weak  [85].  BCPs  indicating  Al-C  bonds  have  low  values  of  pb,  small  negative  values  for 
Hb,  a  kinetic  energy  per  density  ratio  ( Gb/pb )  close  to  one  and  positive  values  for  V2pb. 
These  are  all  indicators  of  bonding  with  metal  [86],  [87].  Additionally,  for  all  but  the  Cp 
variant,  less  than  five  BCP  connecting  each  aluminum  atom  to  its  associated  ligand  were 
found.  This  is  a  key  characterization  of  the  hapticity  of  the  monomer.  Hapticity  is  a 
measure  of  the  position  of  the  aluminum  atom  when  looking  down  at  the  plane  of  the 
ligand.  A  hapticity  of  t]5  would  indicate  shared  bonding  with  all  five  carbon  atoms, 
resulting  in  a  position  of  the  aluminum  atom  in  the  center  of  the  carbon  ring.  At  the  other 
extreme,  a  hapticity  of  77 1  would  indicate  bonding  to  only  one  of  the  carbon  atoms  and 
the  aluminum  atom  would  be  positioned  very  nearly  right  on  top  of  that  atom.  The  fact 
that  QTAIM  analysis  shows  less  than  five  BCPs  between  aluminum  atoms  and  respective 
carbon  rings  would  seem  to  suggest  that  these  monomers  are  not  77s.  As  pointed  out  by 
Farrugia  and  co-workers,  however,  BCPs  are  unreliable  indicators  of  fonnal  hapticity  for 
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metals  bonded  to  7r-carbocyclic  ligands  [88].  A  relatively  high  ellipticity,  e,  around  these 
BCPs  indicates  they  are  in  an  area  of  a  relatively  flat  density,  making  it  difficult  to 
pinpoint  actual  CPs.  This  fact,  as  well  as  the  calculated  geometrical  structures,  indicates 
that  the  formal  hapticity  for  ligands  of  all  variants  is  indeed  rjs,  as  would  be  expected 
from  the  structure  [89].  Average  metrics  for  pertinent  BCPs  for  all  variants  are  listed  in 
Table  8. 


Table  8.  Topological  average  BCP  metrics  for  all  variants. 


Cluster 

BCP  type 

P 

V2p 

£ 

|l\ 

C/P 

H 

AI4CP4 

Al-Al 

0.0378 

-0.0201 

0.4400 

2.7954 

0.1677 

-0.0114 

Al-C 

0.0348 

0.0874 

5.1460 

1.2596 

0.8497 

-0.0077 

ALtCp  4 

Al-Al 

0.0375 

-0.0183 

0.5900 

2.6580 

0.1854 

-0.0115 

Al-C 

0.0393 

0.1193 

12.1500 

1.2228 

0.9708 

-0.0085 

H-H 

0.0041 

0.0131 

0.6739 

0.7330 

0.6240 

0.0007 

H-Cs 

0.0043 

0.0146 

2.1750 

0.7050 

0.6317 

0.0008 

ALjCp  4 

Al-Al 

0.0352 

-0.0154 

0.5033 

2.6057 

0.1804 

-0.0102 

Al-C 

0.0403 

0.1205 

10.0750 

1.2282 

0.9678 

-0.0089 

H-H 

0.0055 

0.0184 

1.8728 

0.7534 

0.6564 

0.0009 

H-C* 

0.0062 

0.0214 

1.3333 

0.7523 

0.6640 

0.0010 

Al4Cp*Ht4 

Al-Al 

0.0353 

-0.0153 

0.4933 

2.5953 

0.1822 

-0.0103 

Al-C 

0.0411 

0.1242 

4.6325 

1.2265 

0.9780 

-0.0091 

H-H 

0.0062 

0.0213 

1.8215 

0.7609 

0.6914 

0.0010 

H-C* 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Al4Cp*lFr4 

Al-Al 

0.0329 

-0.0122 

0.4083 

2.5212 

0.1781 

-0.0089 

Al-C 

0.03895 

0.10833 

15.798 

1.24925 

0.9219 

-0.0090 

H-H 

0.0063 

0.0207 

0.8416 

0.7729 

0.6737 

0.0010 

H-C* 

0.0074 

0.0265 

1.4044 

0.7756 

0.7290 

0.0012 

Al4Cp*Fr4 

Al-Al 

0.0352 

-0.0153 

0.4933 

2.5955 

0.1815 

-0.0102 

Al-C 

0.0410 

0.1237 

4.5800 

1.2275 

0.9756 

-0.0091 

H-H 

0.0062 

0.0211 

1.7895 

0.7623 

0.6913 

0.0010 

H-C* 

0.0040 

0.0140 

1.3475 

0.7070 

0.6868 

0.0008 

*  All  numbers  are  in  atomic  units. 
5  Inter-ligand  BCPs  only. 
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The  information  gained  through  QTAIM  analysis  of  these  clusters  gives  insight 
into  the  bonding  within  the  tetramers.  There  is  a  not  a  simple  trend  as  the  steric  bulk  of 
the  Cp-type  ligand  increases.  It  is  reasonable  to  expect  that  the  addition  of  steric  bulk  to 
the  ligands  of  these  types  of  clusters  would  reduce  tetrameric  stability,  but  the  opposite  is 
the  case.  The  simplest  metric  for  determining  relative  Al-Al  bond  strength  is  the  average 
bond  distance  of  the  six  bonds  that  join  the  core  aluminum  atoms  together.  The  Cp  and 
Cp  variants  have  the  smallest  average  Al-Al  bond  distance  at  2.722  A  and  2.720  A 
respectively,  a  mere  0.19  picometers  apart  and  about  4.5  picometers  from  the  next  closest 
variant.  This  would  seem  to  indicate  that  the  Al-Al  bonds  within  the  core  of  the  Cp  and 
Cp  are  the  strongest  among  all  cluster  variants,  a  notion  in  agreement  with  experimental 
observation  of  only  tetrameric  species  in  solution  for  these  two  compounds.  The  amount 
of  electron  density  at  a  BCP  joining  two  atoms  can  indicate  relative  bond  strength  for 
similar  type  bonds  [85].  Based  on  AIM  analysis,  all  the  Al-Al  bonds  for  all  cluster 
variants  are  shared  or  covalent  type  bonds  with  very  similar  metrics.  As  expected,  the 
total  electron  density  at  the  six  Al-Al  BCPs  is  highest  for  the  Cp  and  Cp  variants.  Table 
9  shows  a  list  of  average  Al-Al  bond  distance  and  total  electron  density  at  the  BCPs 
associated  with  these  bonds  for  all  six  cluster  variants  as  well  as  two  additional 
configurations:  the  Cp  Pl  variant  in  the  configuration  shown  in  Figure  17(b)  and  the  Cp  lPr 
variant  shown  in  the  configuration  of  Figure  17(a).  These  configurations  were  included 
because  of  their  varied  Al-Al  bond  distance. 


Table  9.  Al-Al  bond  distances  and  total  electron  density  at  BCPs  for  eight 

cluster  variants. 


Cp 

O 

T3 

> 

* 

Cp 

Cp*Et 

Cp*lPr 

Cp*Pr 

Cp*Pr  $ 

Cp*ipr  $$ 

Al-Al  dist. 

2.7215 

2.7196 

2.7684 

2.7672 

2.8240 

2.7670 

2.7832 

2.8631 

Pb 

0.2269 

0.2249 

0.2114 

0.2116 

0.1974 

0.2114 

0.2074 

0.1876 

*  Numbers  are  in  A  or  A'3. 

5  Orientation  show  in  Figure  17(b). 
5S  Orientation  show  in  Figure  17(a). 
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There  is  a  clear  correlation  between  the  amount  of  electron  density  at  the  six 
BCPs  joining  the  core  aluminum  atoms  and  the  Al-Al  average  bond  distance.  Figure  24 
shows  a  linear  fit  of  a  plot  of  Al-Al  bond  distance  against  total  electron  density  has  an  R2 
value  of  over  0.99.  This  supports  the  conclusion  that  the  Al-Al  bonds  at  the  core  of  the 
Cp  and  Cp  variants  are  the  strongest.  These  bonds,  however,  are  not  the  only  ones 
involved  in  overall  cluster  stability.  The  H-H  and  H-C  bonds  found  in  all  substituted  Cp 
derivatives  clearly  indicate  there  is  additional  stability  via  non-covalent  ligand 
interactions,  both  in  silico  and  in  solution. 


Figure  24.  Al-Al  bond  distance  vs.  total  Vpb  associated  with  Al-Al  BCPs. 


The  relevant  metric  for  assessing  the  contributions  to  stability  from  the  inter¬ 
ligand  H-H  and  H-C  inter-ligand  bonds  is  the  magnitude  of  electron  density  at  the 
associated  BCPs.  Like  pb  associated  with  Al-Al  bonds,  pb  associated  with  these  closed 
shell  bonds  can  also  indicate  bond  strength  [90],  [91].  It  is  important  to  realize,  however, 
that  similar  amounts  of  electron  density  at  BCPs  associated  with  different  types  of  bonds 
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most  likely  does  not  indicate  that  the  bonds  are  of  the  same  strength  [86].  The  inter¬ 
ligand  H-H  and  H-C  bonds  in  these  clusters  are  closed  shell  bonds  as  opposed  to  the 
shared  Al-Al  bonds  in  the  core.  The  pb  associated  with  these  different  types  of  bonds  will 
have  different  impacts  on  stability.  Additionally,  we  would  not  expect  the  pb  associated 
with  H-H  and  H-C  bonds  to  affect  stability  in  the  same  way.  It  follows  that  there  are  three 
types  of  bonds  that  contribute  to  stability  in  these  types  of  clusters:  Al-Al,  inter-ligand  H- 
H,  and  inter-ligand  H-C  bonds.  The  total  amount  of  electron  density  associated  with  each 
type  of  bond  for  each  cluster  variant  is  listed  in  Table  10.  These  metrics  are  listed  along 
with  the  corresponding  BSSE  corrected  A Htet.  No  solvent  corrections  are  used  so 
fundamental  trends  can  be  analyzed.  Also  included  are  the  same  two  additional 
configurations  examined  in  the  Al-Al  bond  distance  analysis.  The  value  of  A Htet  for 
these  configurations  provides  more  overall  variance  for  subsequent  regression  analysis. 


Table  10.  Total  electron  density  at  Al-Al,  H-H  and  H-C  BCPs  for  various 

clusters. 


Total  Al-Al  pb 

Total  H-H  pb 

Total  H-C  pb 

A Htet  (calc) 

AfiCp4 

0.2269 

0.0 

0.0 

-108.6 

AfiCp  4 

0.2249 

0.0738 

0.0086 

-161.1 

AfiCp  4 

0.2114 

0.1066 

0.0106 

-154.5 

Al4Cp*Et4 

0.2116 

0.1249 

0 

-157.9 

Al4Cp*lFr4 

0.1974 

0.1193 

0.0669 

-162.7 

Al4Cp*Fr4 

0.2114 

0.1242 

0.0159 

-165.7 

(Al4Cp*Fr4)$ 

0.2074 

0.1199 

0.0144 

-160.0 

(AfiCp*1^ 

0.1876 

0.1255 

0.0382 

-135.6 

*  Numbers  are  in  <?-A'3  or  kJ/mol 


$  Orientation  shown  in  Figure  17(b). 
$s  Orientation  shown  in  Figure  17(a). 


Previous  work  by  Boyd  and  co-workers  demonstrated  a  linear  relationship 
between  pb  and  bond  energy  for  various  types  of  bonds  [92].  A  multivariable  linear 
regression  of  the  data  in  Table  10  was  perfonned.  It  was  based  on  the  summed  electron 
density  at  Al-Al,  H-H,  and  H-C  BCPs  to  gain  a  sense  of  their  relative  contributions  to 


tetramer  stability.  The  result  is  this  equation: 
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&Htet  —  a0  +  ax Pai-ai  +  a2  Ph-h  +  a3  Ph-c>  (8) 

where  pX-x  is  the  total  electron  density  at  the  BCPs  associated  with  X-X  bonds  for  a 
given  cluster.  Figure  25  shows  a  plot  of  A Htet  calculated  with  this  equation  against  A Htet 
calculated  through  DFT  vibrational  frequency  analysis.  The  electron  density  associated 
with  H-C  inter-ligand  bonds,  governed  by  the  coefficient  a3,  contributes  the  least  to 
tetramer  stability.  The  optimal  fit  gives  ratios  %/ a3  ~  3.8  and  a2/a3  ~  1.3.  Thus, 
approximately  65%  of  the  variation  in  tetramer  binding  energy  can  be  explained  by  the 
electron  density  at  Al-Al  BCPs,  which  varies  linearly  with  Al-Al  separation.  The 
remainder  can  be  attributed  to  inter-ligand  effects. 


Figure  25.  AHtet  calculated  via  DFT  frequency  analysis  compared  to  A Htet 

calculated  via  equation  (8). 
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This  analysis  helps  explain  why  the  Cp  variant  has  the  lowest  binding  energy  of 
all  clusters  considered  here,  despite  having  the  second  smallest  Al-Al  bond  distance  and 
the  largest  pb  associated  with  Al-Al  bonds.  The  Cp  variant  has  no  inter-ligand  H-H  or  H- 
C  bonds  and  therefore  no  additional  contribution  to  tetramer  stability.  Compared  to  the 
Cp  variant,  the  Cp  variant  has  only  a  slightly  smaller  average  Al-Al  bond  distance  and 
slightly  less  pb  associated  with  Al-Al  bonds,  but  because  of  the  additional  pb  associated 
with  inter-ligand  H-H  and  H-C  bonds,  the  Cp  variant  is  52.5  kJ/mol  more  stable.  This  is 
also  why  the  Cp  Pr  variant  is  slightly  more  stable  than  the  Cp*  despite  having  the  same 
amount  of  pb  associated  with  Al-Al  bonds.  The  Cp*Pl  has  roughly  20%  more  electron 
density  associated  with  inter-ligand  H-H  and  H-C  bonds.  It  is  clear  that  inter-ligand  H-H 
and  H-C  bonds  make  non-trivial  contributions  to  the  stability.  More  steric  bulk  does  not 
necessarily  mean  more  stability,  however.  Even  though  the  Cp  variant  has  more  steric 
bulk  than  the  Cp  variant,  the  Cp  variant  is  slightly  more  stable.  Despite  the  addition  of 
four  methyl  groups  to  each  ligand  in  a  Cp  variant,  the  average  Al-Al  distance  in  the 
resulting  Cp  variant  is  virtually  identical  (within  0.2  picometers)  to  that  of  the  Cp 
variant.  The  addition  of  a  fifth  methyl  group,  however,  results  in  an  average  Al-Al 
distance  in  a  Cp*  variant  that  is  5.5  picometers  longer  than  that  of  a  Cp  variant.  While  the 
Cp*  variant  has  more  pb  associated  with  inter-ligand  bonds  than  the  Cp  variant,  because 
of  its  extra  steric  bulk,  the  Cp  variant  has  more  pb  associated  with  Al-Al  bonds  due  to  its 
smaller  average  Al-Al  distance.  The  extra  pb  associated  with  inter-ligand  bonds  in  the 
Cp*  variant  is  not  sufficient  to  compensate  for  the  loss  of  pb  associated  with  Al-Al  bonds. 
Ultimately  there  is  no  simple  correlation  between  ligand  bulk  and  enthalpic  stability. 

D.  EQUILIBRIUM  CONSTANT  CALCULATIONS 

To  further  investigate  the  behavior  of  this  series  of  aluminum  tetramer  clusters, 

their  monomer/tetramer  equilibrium  in  solution  was  examined.  Availability  of  free 

monomers  in  solution  is  likely  a  key  element  in  the  fonnation  of  larger  metalloid  clusters 

such  as  Al50Cp*i2.  As  tetrameric  forms  disproportionate  into  trivalent  species  (e.g., 

AlCp  3)  and  bulk  aluminum,  the  presence  of  free  monomers  in  solution  could  allow  for 

their  binding  to  small  particles  of  bulk  aluminum  to  fonn  metalloid  clusters.  To 

understand  the  behavior  of  the  clusters  in  this  work,  analysis  similar  to  that  of  Huber  and 
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Schnockel  following  synthesis  of  the  Cp  variant  was  performed  [38],  The 
thermodynamic  results  for  all  six  variants  were  used  to  calculate  the  equilibrium  constant, 
K ,  at  room  temperature  using  the  van  ‘t  Hoff  equation. 

AH-TAS 

Keq  =  e  .  (9) 

Despite  excellent  agreement  with  tetramerization  enthalpies,  K  values  calculated 
via  DFT  differ  dramatically  from  experiment  due  to  the  exponential  dependence  on  AS. 
As  discussed  previously,  no  computational  method  examined  was  found  to  predict 
entropy  changes  with  sufficient  accuracy  to  match  experimental  K  values.  Schnockel  and 
Huber  linearly  scaled  calculated  A Htet  and  A Stet  values  such  that  the  results  for  the  Cp 
variant  would  match  experimental  values.  Because  of  the  poor  perfonnance  by  DFT  to 
predict  A Stet,  however,  rather  than  scaling  the  values  for  A Htet  and  A Stet  separately,  the 
DFT  results  for  both  A Htet  and  A Stet  are  used  to  calculate  K  for  each  cluster  before  any 
scaling.  The  resulting  K  values  are  then  normalized  so  that  the  result  for  the  Cp  variant 
corresponds  to  a  monomer  to  tetramer  ratio  of  1:449,  the  experimental  ratio.  Ratios  for 
the  other  clusters  can  then  be  calculated  from  their  respective  nonnalized  K  values.  The 
larger  the  K  value,  the  larger  the  ratio  of  tetramers  to  monomers  and  hence  the  fewer  free 
monomers  in  solution.  The  results  of  this  analysis,  as  well  as  experimental  values,  are 
shown  in  Table  11.  It  should  be  stressed  that  this  method  will  likely  not  indicate  which 
clusters  are  likely  to  have  monomers  found  in  solution,  but  instead  gives  a  comparison  of 
which  clusters  are  more  or  less  likely  to  have  monomers  in  solution  as  compared  to  the 
Cp  variant.  Rather  than  attempting  to  produce  predictive  values  for  K,  this  method  is 
comparative. 
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Table  1 1 .  Experimental  and  calculated  equilibrium  constants  at  298. 1 5K. 


AI4CP  4 

AI4CP4 

AfiCp  4 

Al4Cp*Et4 

Al4Cp*Pr4 

Al4Cp*lPr4 

Experimental  Values 

A  HtPt 

-150 

- 

- 

- 

-160 

-158 

AStet 

-0.300 

- 

- 

- 

-0.398 

-0.477 

K 

4.1  x  101U 

- 

- 

- 

1.8  x  107 

5.8  x  102 

Ratio  (mon/tet) 

1:449 

- 

- 

- 

1:65 

1:5 

Calculated  Values 

A  Htet 

-148.4 

-106.7 

-153.1 

-151.5 

-158.7 

-155.8 

AStet 

-0.639 

-0.380 

-0.437 

-0.639 

-0.664 

-0.663 

K 

4.2  x  10'8 

7.0  x  10'2 

9.9  x  10J 

1.5  x  Hf7 

1.3  x  Hf7 

4.6  x  10~* 

Normalized  K  Values 

K 

4.1  x  101U 

6.8  x  1016 

9.6  x  1021 

1.4  x  1011 

1.3  x  1011 

4.5  x  10IU 

Ratio  (mon/tet) 

1:449 

1:16,132 

1:313,091 

1:613 

1:598 

1:460 

*  A Htet  numbers  are  in  kJ/mol;  AStet  numbers  are  in  kJ/mol-K 


The  results  show  that  the  K  values  for  Cp  ,  Cp*lPr,  Cp  Pr  and  Cp*Et  variants  are  all 
within  an  order  of  magnitude  and  should  have  similar  monomer/tetramer  ratios.  The 
results  also  show  that  K  values  of  the  Cp  and  Cp  variants  are  respectively  6  and  1 1 
orders  of  magnitude  larger  than  that  of  the  Cp  variant.  This  suggests  that  there  are  likely 
to  have  significantly  lower  concentrations  of  Cp  and  Cp  monomers  as  compared  to  the 
Cp  variant.  It  is  likely  that  only  differences  of  more  than  an  order  of  magnitude  are 
significant.  It  is  unreasonable  to  differentiate  the  concentrations  of  the  Cp*,  Cp  lPr,  Cp*Pr 
and  Cp*Et  variants  since  calculated  K  values  for  these  variants  are  all  relatively  close. 
Since  monomers  of  the  Cp*  are  known  to  exist  in  solution,  then  based  on  this  analysis,  it 
is  reasonable  to  expect  to  find  monomers  of  the  Cp  lPr,  Cp*Pl  and  Cp  Et  variants  as  well. 
While  the  Cp*Et  variant  has  not  yet  been  synthesized,  monomers  of  both  the  Cp*lPr  and 
Cp  Pr  variants  were,  in  fact,  found  in  solution.  This  analysis  also  indicates  that  monomers 
of  Cp  variant  should  be  somewhat  less  prevalent  in  solution  as  compared  to  the  Cp 
variant.  Monomers  of  the  Cp  variant  were  never  found  in  solution,  though  it  is  important 
to  note  that  disproportiation  of  the  Cp  variant  into  trivalent  species  and  bulk  aluminum 
occurs  significantly  earlier  than  the  other  variants,  at  85°C  below  room  temperature  [38]. 
Finally,  the  analysis  suggests  that  monomers  of  the  Cp  variant  are  of  significantly  lower 
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concentration  than  those  of  the  Cp  variant.  This  is  consistent  with  the  experimental  fact 
that  monomers  of  the  Cp  variant  were  never  found  in  solution  [38]. 

The  behavior  of  these  systems  with  regard  to  monomer/te trainer  equilibrium  may 
be  of  importance  for  alternate  synthesis  methods  such  as  those  reported  by  Roesky  [33] 
and  Fisher  [63],  which  can  produce  AIR  monomers  without  the  need  for  a  complex  co¬ 
condensation  apparatus.  A  bias  towards  monomers  in  solution  at  higher  temperatures  may 
allow  for  greater  synthetic  availability  of  monomers  in  these  solution-based  synthesis 
methods.  However,  it  appears  that  the  complex  balance  between  monovalent  Al-Al 
bonding  in  the  core  and  non-covalent  ligand  interactions  in  the  exterior  do  not  lead  to  a 
general  trend  with  steric  bulk,  but  rather  requires  analysis  of  each  individual  system. 
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III.  LARGER  LIGAND  VARIANTS 


A.  EFFECTS  OF  ADDITIONAL  STERIC  BULK 

While  the  addition  of  steric  bulk  to  the  unsubstituted  Cp  variant  does  not  have  a 
direct  correlation  with  enthalpic  stability,  it  does  appear  that  the  addition  of  steric  bulk 
leads  to  smaller  equilibrium  constants,  and  therefore  larger  monomer/tetramer  ratios  in 
solution.  The  two  clusters  investigated  up  to  this  point  that  have  the  highest  K  also  have 
the  least  amount  of  steric  bulk.  The  Cp  variant,  however,  has  only  slightly  less  steric 
bulk  than  the  Cp  variant,  yet  has  a  significantly  higher  K.  To  further  understand  the 
relationship  between  steric  bulk  and  monomers  in  solution  an  examination  of  four 
additional  Cp-type  variants  is  conducted. 

1.  Al4(C5Me4Ph)4 

In  2005  Buchin  et  al.  reported  possible  synthesis  of  the  cluster  Al4Cp*Ph4  (Cp*Ph  = 
C5Me4Ph),  but  the  cluster  failed  to  crystalize  and  experimental  27A1  NMR  results  were 
inconclusive  [25].  Buchin  et  al.  reported  a  27A1  NMR  chemical  shift  of  102.5  ppm,  which 
is  significantly  higher  than  that  of  the  Cp  and  Cp  tetramers,  and  no  DFT  calculations 
were  reported  [25].  Of  significance  was  the  fact  that  only  one  NMR  signal  was  detected. 
Variable  temperature  NMR  by  David  Morris  at  NAWC  China  Lake  also  confirmed  the 
presence  of  only  one  signal  up  to  a  temperature  of  100  degrees  C.  The  Cp*Ph  variant, 
shown  in  Figure  26,  has  the  most  steric  bulk  of  any  Cp*  derivative  investigated  so  far. 
Assuming  the  single  NMR  signal  is  the  tetramer,  then  the  absence  of  monomers  in 
solution  is  somewhat  surprising.  To  examine  this  cluster,  first  the  structure  of  the  variant 
must  be  determined.  DFT  NMR  analysis  is  then  conducted  and  compared  with 
experimental  results  from  2005  as  well  as  more  recent  results  from  NAWC  China  Lake. 
DFT  vibrational  frequency  analysis  is  then  conducted  to  detennine  thermodynamic 
characteristics  which  are  not  available  experimentally.  Finally,  calculated  thermodynamic 
values  can  be  used  to  detennine  the  value  of  the  equilibrium  constant  for  this  cluster  and 
compared  with  previous  six  clusters  studied. 
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(a) 


Hydrogens  of  tetramer  removed  for  clarity. 


Figure  26.  Monomer  (a)  and  tetramer  (b)  of  Cp  Ph4  variant 


*Ph 


a.  Structure  of  Al4Cp  4 

Following  synthesis  of  the  CpPh  tetramer,  all  attempts  to  crystalize  a  structure 
failed  [25].  No  experimental  structural  data  is  available.  The  ortho  C-H  of  the  phenyl  ring 
protrudes  in  towards  the  aluminum  atom  in  the  energy-minimized  structure.  This  is  in  a 
way  similar  to  the  Cp*lPr  variant,  so  a  tetramer  configured  as  in  Figure  18  was  considered 
to  be  a  likely  global  minimum.  Because  of  the  possibility  of  twisting  of  the  phenyl  ring, 
however,  it  is  not  obvious  that  this  configuration  would  be  the  lowest  energy  structure.  In 
order  to  examine  a  large  number  of  possible  structures,  a  genetic  algorithm  (GA)  was 
implemented  using  a  more  efficient  DFT  code,  GPAW.  GPAW  is  a  DFT  python  code 
based  on  the  projector  augmented  wave  (PAW)  method  [93],  [94].  The  code  operates 
within  the  Atomic  Simulation  Environment  (ASE)  which  provides  a  means  for  creating 
and  manipulating  atomic  structures  [95].  The  PAW,  as  implemented  within  GPAW,  takes 
advantage  of  the  frozen  core  approximation  [96].  This  approximation  leverages  the 
known  behavior  of  electrons  that  are  not  part  of  the  valence.  Core  electrons  are  treated 
differently  and  with  much  less  computational  effort  than  the  valence  electrons.  The 
resulting  core  wavefunctions  are  smoothly  combined  with  the  valence  pseudo 
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wavefunctions  to  arrive  at  an  approximation  of  the  all  electron  Kohn-Sham  wavefunction 
[93].  This  provides  significant  computational  efficiency  over  methods  that  solve  the  KS 
orbital  eigenvalue  problems  for  all  elections,  the  approach  used  by  Gaussian. 
Additionally,  within  the  PAW  framework,  LCAO  can  be  used  to  realize  greater 
computational  efficiency.  Rather  than  using  large  collections  of  plane  waves,  which  are 
typically  used  in  the  PAW  framework,  LCAO  uses  smaller  numbers  of  atomic-like  basis 
functions  to  approximate  wavefunctions  [94].  This  allows  for  faster  calculations  with 
only  a  small  loss  to  overall  accuracy.  Ultimately,  energy  calculations  using  GPAW  will 
most  likely  be  less  accurate  than  those  performed  with  Gaussian,  but  they  will  be  done  in 
far  less  time.  Additionally,  the  errors  suffered  within  GPAW  are  systematic,  meaning 
they  would  be  common  amongst  similar  structures  and  especially  among  like  tetramers 
with  monomers  in  different  configurations.  For  this  reason,  GPAW  was  used  to 
implement  the  GA  designed  to  find  the  lowest  energy  tetramer.  To  ensure  accurate 
dispersion  effects  were  captured,  a  critical  component  of  these  clusters,  the  dispersion 
correction  method  of  Tkatchenko  and  Scheffler,  as  implemented  within  GPAW,  was 
used  [97]. 

To  create  a  tetramer,  a  single  monomer  is  loaded  into  the  ASE  environment  and 
copied  three  times.  The  three  copies  are  rotated  and  translated  in  such  a  way  as  to 
produce  an  ideal  tetrameric  aluminum  core  with  a  given  Al-Al  interatomic  distance.  For 
the  case  of  the  Cp  Ph  variant  the  chosen  distance  was  2.8  A,  roughly  the  average  Al-Al 
distance  for  the  Cp  lPr  variant.  All  four  monomers  are  then  individually  rotated  about  the 
axis  through  their  respective  Cp  ring  center.  Each  is  rotated  a  random  number  of  degrees 
and  the  number  of  degrees  rotated  is  recorded  into  bins  of  10  degrees.  A  set  of  four  bins 
would  describe  the  configuration  of  the  tetramer.  For  example,  a  tetramer  with  a 
configuration  of  [0,  2,  18,  35]  would  represent  a  tetramer  whose  four  monomers  were 
rotated  between  0-10  degrees,  20-29  degrees,  180-189  degrees  and  350-359  degrees 
respectively.  Order  is  important.  The  configuration  [0,  18,  35,  2]  would  be  an  entirely 
different  structure  due  to  the  order  of  rotation  operations  and  the  interaction  of  each 
monomer  with  its  neighbors.  The  geometry  of  the  created  structure  is  then  optimized 
using  GPAW.  This  optimization  tends  to  keep  the  monomers  in  their  initial  randomly 
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selected  configuration,  rotating  only  slightly  before  a  local  potential  energy  minimum  is 
reached.  Attempting  to  systematically  check  all  possible  monomer  configurations  would 
require  optimizing  1,679,616  structures.  The  optimization  of  100  structures  takes  about 
1,800  hours  with  72  CPUs.  On  a  super  computing  cluster  with  multiple  nodes  and 
thousands  of  CPUs  these  structures  can  be  optimized  simultaneously.  However,  the  cost 
of  optimizing  all  100  structures  would  be  about  130,000  CPU  hours,  which  represents 
more  than  10%  of  the  typical  annual  allotment  of  CPU  time  for  research  groups  using 
DOD  HPC  resources.  Investigating  just  1%  of  the  total  possible  structures  would  require 
nearly  22  million  hours  of  CPU  time  and  more  than  likely  would  only  give  about  a  1% 
chance  of  finding  the  global  minimum.  In  order  to  find  the  global  minimum  for  the  Cp*Ph 
tetramer  in  a  practical  amount  of  CPU  time,  a  GA  in  the  spirit  of  Deaven  and  Ho  was 
used  [98].  There  are  six  basic  steps  for  this  GA.  First,  create  an  initial  group  of  randomly 
generated  candidates.  Second,  select  from  these  candidates  the  lowest  energy  structures 
to  form  an  initial  population.  Third,  preferentially  select  candidates  from  the  population 
for  pairing  and  create  a  new  offspring.  Fourth,  mutate  the  offspring  with  some  non-zero 
probability.  Fifth,  relax  the  offspring  and  include  in  the  current  population  if  the  resulting 
energy  is  low  enough.  Finally,  repeat  the  process  until  the  global  minimum  has  likely 
been  found.  The  application  of  this  process  for  the  Cp  Ph  variant  borrows  heavily  from  the 
work  of  Vilhelmsen  and  Hammer  on  metal  clusters  [99], 

(1)  Starting  Candidates 

The  initial  set  of  randomly  generated  candidates  should  be  as  diverse  as  possible, 
but  there  is  no  rule  of  thumb  as  to  how  large  the  list  should  be.  Vilhelmsen  and  Hammer 
were  studying  groups  of  metal  clusters  with  less  than  50  atoms.  Random  positions  for 
each  atom  were  chosen  to  create  initial  structures.  This  results  in  a  significant  number  of 
degrees  of  freedom  for  each  cluster,  however  it  is  very  likely  that  several  different 
configurations  could  be  optimized  into  the  global  minimum  as  there  are  no  constraints 
holding  an  atom  in  place  other  than  the  existence  of  neighboring  atoms.  The  creation  of  a 
Cp  Ph  cluster  only  has  four  degrees  of  freedom,  but  it  is  much  less  likely  for  many 
different  configurations  to  relax  into  the  global  minimum.  For  this  reason,  rather  than  an 

initial  set  of  20  candidates,  as  recommended  by  Vilhelmsen  and  Hammer,  a  set  of  50  was 
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chosen  [99].  To  prevent  creating  similar  candidates,  the  configuration  of  each  structure, 
in  terms  of  angular  bins,  was  recorded  and  repeat  configurations  were  rejected. 

(2)  Creating  Initial  Population 

A  small  population  allows  stronger  candidates  to  mate  more  often,  but  lacks 
diversity.  Vilhelmsen  and  Hammer  found  that  small  populations  required  less  generations 
to  reach  a  minimum,  but  the  resulting  minimum  was  less  likely  to  be  the  actual  global 
minimum  [99].  A  small  population  size  was  considered  10  candidates,  while  a  large 
population  was  40  candidates  [99].  For  this  GA,  an  initial  population  of  15  was  chosen, 
as  speed  was  of  primary  importance.  Additionally,  it  is  not  necessary  to  find  the  actual 
global  minimum.  As  long  as  the  resulting  candidate  was  within  a  few  kilojoules  per  mole 
of  the  actual  global  minimum  then  this  would  represent  sufficient  accuracy  for  ultimate 
monomer/te trainer  equilibrium  analysis.  A  population  of  the  15  lowest  energy  candidates 
among  the  entire  set  of  structures  is  maintained  throughout  the  entire  GA  run. 


(3)  Selection  of  Mates 


Within  the  population,  it  is  important  to  preferentially  select  the  most  fitting 
candidates  for  mating.  Randomly  choosing  two  candidates  will  ultimately  slow  the  GA’s 
progression  towards  a  global  minimum.  It  is  important  to  select  a  candidate  based  not  on 
position  in  the  population,  but  on  actual  strength  of  energy  compared  to  the  rest  of  the 
population.  For  this  reason,  the  fitness  equation  used  by  Vilhelmsen  and  Hammer  was 
adopted  [99]. 

Fi  =  1~[  1  -  tanh(2p,  -  1)],  with  Pi  =  Ei~EmJn  .  (10) 

z  ^ max  min 

Here  Emin  is  the  energy  of  the  lowest  energy  candidate  in  the  population  and  Emax  is  the 
energy  of  the  candidate  with  the  highest  energy  in  the  population.  The  resulting  fitness, 
Fj,  ranges  from  0.88  to  0.12  and  is  further  adjusted  by  the  number  of  times  a  candidate 
has  been  selected  for  mating,  nt  [99].  The  final  probability  of  selection  is: 


P  =  F  — - — 

1  1  -JE+ni 


(11) 


This  construct  preferentially  selects  the  lowest  energy  candidates,  but  ensures  that 
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particular  candidates  are  not  picked  too  many  times.  This  helps  to  maintain  diversity  in 
the  population. 


(4)  Pairing 

Pairing  of  two  selected  mates  is  straightforward  considering  there  are  only  four 
degrees  of  freedom  for  each  structure.  Each  mate  must  contribute  at  least  one  monomer 
to  the  offspring.  This  results  in  only  14  possible  combinations  of  monomers  for  the 
offspring.  During  a  pairing,  one  of  the  14  possible  combinations  is  randomly  selected  and 
a  new  structure  is  created.  The  configuration  of  the  structure  is  derived  from  the 
appropriate  configuration  numbers  of  the  two  mates.  For  example,  two  mates  with 
configurations  [0,  5,  18,  30]  and  [1,  6,  19,  31]  could  create  offspring  with  configurations 
of  [0,  5,  19,  31]  or  [1,  5,  19,  31],  but  not  one  with  a  configuration  of  [18,  19,  30,  31].  The 
resulting  structure  cannot  be  accepted  until  its  configuration  is  compared  to  the 
configurations  of  all  structures  in  the  entire  GA.  If  the  configuration  is  a  repeat  the 
structure  is  rejected  and  another  one  of  the  14  possible  combinations  for  the  two  selected 
mates  is  tried.  In  the  case  of  a  mutation,  the  selected  configuration  is  not  eliminated  and 
can  be  tried  again.  This  continues  until  all  14  combinations  are  exhausted  at  which  point 
the  two  mates  are  considered  incompatible  and  new  mates  are  chosen. 

(5)  Mutation 

Mutations  are  a  critical  aspect  of  any  GA.  They  introduce  random  diversity  into 
the  existing  population.  With  no  mutations,  there  is  a  risk  of  the  entire  population  falling 
into  a  wide  local  minimum  resulting  in  the  GA’s  inability  to  find  the  global  minimum. 
Mutating  too  often,  however,  introduces  too  much  randomness  into  the  population  and 
adversely  affects  the  GA’s  ability  to  progress  towards  the  global  minimum.  For  this 
work,  a  30%  mutation  rate  is  used  to  match  the  rate  used  by  Vilhelmsen  and 
Hammer  [99].  To  accomplish  this  each  of  the  four  monomers  are  tested  for  mutation  in 
turn.  There  is  an  8.5%  chance  that  a  monomer  would  mutate,  which  results  in  an  overall 
chance  of  any  one  of  the  four  monomers  mutating  of  30%.  Additionally,  it  is  possible  for 
more  than  one  monomer  to  mutate.  The  mutation  is  a  random  rotation  about  the  axis 
through  the  Cp  center  of  the  mutated  monomer.  The  configuration  of  the  new  structure  is 
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then  updated  based  on  which  angular  bin  the  mutated  monomer  would  ultimately  fall  in 
after  rotation. 

(6)  Optimization  and  Entry  into  Population 

All  structures  are  optimized  to  the  nearest  local  minimum  using  GPAW.  In  order 
to  take  advantage  of  the  volume  of  CPUs  available  on  DOD  HPC  clusters,  it  is  desirable 
to  optimize  as  many  structures  simultaneously  as  possible.  However,  in  an  effort  to 
extract  the  maximum  efficiency  of  the  GA,  a  departure  from  the  approach  of  Vilhelmsen 
and  Hammer  is  necessary  to  ensure  a  maximum  number  of  generations.  Therefore,  a 
generational  approach  similar  to  that  of  Johnston’s  Birmingham  cluster  GA  is 
employed  [100].  A  small  group  of  offspring  are  selected  and  allowed  to  fully  optimize.  A 
new  population  is  created,  including  newly  relaxed  offspring  if  appropriate.  The  members 
of  this  new  population  are  then  given  fitness  scores  and  another  generation  is  selected. 
While  there  is  no  rule  of  thumb  as  to  the  number  and  size  of  generations  necessary  for 
optimal  efficiency,  a  generation  size  of  ten  was  chosen  to  allow  for  ten  total  generations 
by  the  time  150  structures  were  relaxed.  This  does  not  include  generation  0,  which  would 
be  the  population  created  from  the  first  fifty  randomly  created  structures. 

(7)  Tennination  of  the  GA  and  Optimization  in  Gaussian 

The  goal  was  for  the  GA  to  narrow  in  on  the  global  minimum  within  the  first  ten 
generations,  or  150  total  structures  relaxed.  For  the  CpPh  variant,  the  lowest  energy 
structure  was  produced  during  the  fifth  generation.  By  the  end  of  the  tenth  generation  the 
top  ten  structures,  though  initially  configured  differently,  all  optimized  to  virtually  the 
same  configuration.  This  configuration  is  shown  in  Figure  17(c)  with  Ph  groups  in  place 
of  the  Pr  groups.  This  is  not  surprising  as  it  is  logical  that  the  Ph  groups  would  prefer  to 
settle  into  the  gaps  between  adjacent  ligands  in  a  manner  similar  to  that  of  the  iPr  groups 
in  the  CplPr  variant.  While  it  is  impossible  to  know  if  the  actual  global  minimum  was 
found  without  an  experimental  crystal  structure,  based  on  the  inclusion  of  mutations  to 
prevent  stagnation  it  is  likely  a  structure  very  close  to  the  global  minimum,  if  not  the 
global  minimum,  was  in  fact  found.  It  is  also  clear  that  there  are  many  different  initial 
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configurations  that  would  be  capable  of  relaxing  into  the  global  minimum,  making  it 
easier  to  find. 

To  examine  the  performance  of  the  GA,  all  of  the  randomly  generated  structures 
produced  for  the  Cp*Ph  GA  run  as  well  as  the  GA  runs  for  the  three  iPr  variants  discussed 
in  the  next  section  were  compared  to  offspring  produced  by  the  GA.  Histograms  of  both 
the  random  sampling  and  GA  offspring  were  created  and  are  shown  in  Figures  27  and  28 
respectively.  Relative  energies  with  respect  to  the  lowest  energy  structure  found  for  each 
respective  variant  are  plotted.  The  random  sampling  has  a  roughly  nonnal  distribution,  as 
would  be  expected,  and  also  fails  to  find  the  global  minimum  for  any  variant.  The 
histogram  of  the  GA  offspring  shows  the  remarkable  ability  of  this  approach  to  find 
significant  numbers  of  low  energy  structures,  increasing  the  chances  of  finding  the  global 
minimum.  It  is  also  important  to  note  that  more  than  95  percent  of  the  randomly 
generated  structures  were  more  than  25  kJ/mol  higher  in  energy  then  the  minimum  found 
by  the  GA.  Conversely,  nearly  60  percent  of  the  structures  created  by  the  GA  were  within 
25  kJ/mol  of  the  eventual  minimum  found.  Ideally,  the  GA  approach  would  be  validated 
against  a  variant  with  an  experimentally  known  structure.  Unfortunately,  there  are  only 
two  variants  for  which  there  is  experimental  structural  data,  the  Cp  and  the  Cp  variants. 
The  Cp  is  symmetric  about  rotation,  so  there  would  be  thousands  of  degenerate 
configurations  which  would  likely  be  found  quite  easily  by  just  a  random  sampling  of 
structures.  Similarly,  as  previously  mentioned,  the  Cp  is  likewise  insensitive  to  monomer 
rotation  and  therefore  would  not  be  a  good  test  for  the  GA. 
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Figure  27.  Histogram  of  calculated  energy  for  randomly  created  structures. 


Figure  28.  Histogram  of  calculated  energy  for  structures  produced  by  GA. 


A  sampling  of  structures  from  the  GA  were  optimized  in  Gaussian  with  M06-2X 
and  the  split-valence  double  zeta  basis  set  6-31G(d).  Results  showed  that  GA  structures 
within  4-5  kJ/mol  of  each  other  were  essentially  degenerate  when  optimized  with 
Gaussian,  even  if  they  did  not  have  the  same  configuration.  Relative  ranking  for 
structures  with  larger  energy  differences,  however,  were  preserved.  Further  optimization 
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of  the  lowest  energy  structure  for  the  Cp  Ph  variant  was  done  using  the  triple  zeta  cc- 
pVTZ  basis  set. 

b.  NMR  Results 

The  inability  to  obtain  crystals  of  the  Cp  Ph  tetramer  forces  the  reliance  on  27A1 
NMR  to  identify  the  compound  in  solution.  The  Cp*  and  Cp  tetramers  have  27A1  NMR 
chemical  shifts  of  -82  ppm  and  -8 1  ppm  respectively.  They  are  the  only  two  tetramers  for 
which  crystals  have  been  obtained.  The  Cp*Pr  and  Cp*lPr  variants  never  crystalized,  but 
compelling  thermochemical  and  DFT  data  suggest  that  the  experimental  signals  of  -86 
ppm  and  -89  ppm  were  in  fact  the  respective  tetramers  of  these  two  variants.  This  would 
seem  to  indicate  is  that  the  27A1  NMR  chemical  shift  for  the  Cp*Ph  variant  should  be 
relatively  close  to  these  numbers.  The  shift  of  -102.5  found  by  Buchin  et  al.  in  2005  was 
too  far  from  the  signals  of  the  known  Cp  and  Cp  tetramers  for  them  to  definitively  say 
that  they  had  synthesized  the  Cp*Ph  tetramer  [25].  David  Morris  from  NAWC,  China 
Lake,  attempted  to  synthesize  the  Cp  Ph  tetramer  using  multiple  routes  and  consistently 
detected  a  signal  at  -82  ppm.  Buchin  et  al.  used  a  synthesis  process  similar  to  that  of 
Roesky  [33],  but  with  the  halogen  bromine  instead  of  chlorine  to  improve  yield  [25]. 
Morris  used  the  same  synthesis  method  as  Buchin  et  al.  as  well  as  the  method  of  Fischer 
[63]  reported  on  in  2013.  It  is  unclear  as  to  the  discrepancy  in  27NMR  results,  but  DFT 
calculated  NMR  chemical  shifts  predict  a  shift  of  -84.7  ppm,  indicating  Morris  has  more 
than  likely  synthesized  the  tetramer.  Additionally,  variable  temperature  NMR  data  up  to 
100  degrees  C  contains  only  one  signal  at  -82  ppm,  indicating  no  monomers  in  solution. 

c.  Thermochemistry  and  Equilibrium  Constant 

Vibrational  frequency  analysis  of  the  Cp*Ph  variant  indicates  that  it  is  more 
enthalpically  stable  than  the  other  six  variants  previously  investigated.  Calculated  A Htet 
for  the  Cp*Ph  variant  was  -182.9  kJ/mol.  This  is  significantly  lower  than  the  A Htet  of  the 
Cp  Pr  variant,  -158.7  kJ/mol.  The  lower  value  of  A Htet  for  the  Cp*Ph  variant  is  consistent 
with  its  average  Al-Al  interatomic  distance  and  with  QTAIM  analysis.  The  Cp*Ph  has  a 
longer  average  Al-Al  bond  distance,  and  less  pb  associated  with  Al-Al  BCPs  as 
compared  with  the  Cp*Pr  variant.  Total  pb  associated  with  inter-ligand  H-C  and  H-H 
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NCIs,  however,  is  0.1640  A"3  for  the  Cp  Ph  variant.  This  is  almost  20  percent  more  than 
the  total  inter-ligand  pb  calculated  for  the  Cp1’1  variant.  Table  12  shows  these  metrics  for 
both  variants.  The  Cp  Ph  variant  is  another  demonstration  that  weakening  of  the  bonding 
in  the  aluminum  core  of  a  tetramer  due  to  ligand  hindrance  can  be  offset  by  additional 
inter-ligand  NCIs. 


Table  12.  Average  Al-Al  distance  and  pb  for  the  Cp  Pl  and  Cp*Ph  variants. 


Al-Al  dist 

Total  Al-Al  pb 

Total  H-H  pb 

Total  H-C  pb 

Al4Cp*Pr4 

2.7670 

0.2114 

0.1242 

0.0159 

Al4Cp*Ph4 

2.8052 

0.2046 

0.0934 

0.0706 

*  Numbers  are  in  A  or  A"3. 


Available  thermochemical  data  now  allows  for  an  assessment  of  the  relative 
monomer  to  tetramer  ratio  in  solution  as  compared  to  the  other  six  variants.  Figure  29 
shows  a  logarithmic  plot  of  the  monomer  to  tetramer  ratio  of  the  six  variants  previously 
discussed  as  well  as  the  Cp  Ph  variant.  All  numbers  were  normalized  to  produce  a 
monomer  to  tetramer  ratio  of  1 :449  for  the  Cp  variant,  the  known  experimental  value. 
The  results  show  that  the  Cp  Ph  variant  is  distinct  from  the  Cp*,  Cp  Et,  Cp  Pl  and  Cp  lPr 
variants,  which  all  have  roughly  the  same  monomer  to  tetramer  ratio.  The  Cp*Ph  variant  is 
also  distinct  from  the  Cp  variant  which  has  a  much  smaller  monomer  to  tetramer  ratio. 
The  Cp*Ph  has  a  monomer  to  tetramer  ratio  almost  identical  to  that  of  the  Cp  variant,  a 
variant  for  which  monomers  are  not  found  in  solution.  Due  to  reasons  previously 
discussed,  this  type  of  analysis  can  only  inform  of  relative  trends.  The  analysis  suggests 
that  the  Cp*Ph  variant  is  likely  to  have  a  significantly  smaller  monomer  to  tetramer  ratio 
compared  to  the  Cp  ,  Cp  Et,  Cp  Pl  and  Cp*lPr  variants  and  a  significantly  higher  monomer 
to  tetramer  ratio  as  compared  with  the  Cp  variant.  Based  on  experimental  equilibrium 
data  for  the  Cp*,  Cp*Pl,  Cp  lPr,  Cp  and  CpA  variants,  a  definitive  conclusion  with  respect  to 
monomers  in  solution  for  the  Cp*Ph  variant  cannot  be  made.  However,  the  results 
calculated  with  DFT  and  shown  in  Figure  29  are  consistent  with  the  fact  that  monomers 
of  the  Cp*Ph  were  not  found  in  solution  by  two  different  experimental  groups,  even  at 
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elevated  temperatures.  The  fact  that  the  calculated  ratio  for  the  CpPh  variant  would  fall  so 
much  lower  than  the  ratio  for  the  Cp*,  Cp  Et,  Cp  Pl  and  Cp  lPr  variants  is  surprising 
considering  that  the  Cp  Ph  carries  with  it  significantly  more  bulk.  As  QTAIM  analysis  has 
shown,  though,  this  bulk  is  not  steric  and  makes  non-trivial  contributions  to  stability  that 
bind  the  Cp  Ph  variant  more  strongly  than  other  less  bulky  variants.  Adding  bulk  is  not 
sufficient  to  increase  the  monomer/tetramer  ratio  in  solution  due  primarily  to  the  small 
effect  on  Al-Al  distance.  It  is  likely  that  continued  stretching  of  the  Al-Al  distance  with 
more  ligand  hindrance  is  needed  to  degrade  aluminum  bonding  in  the  core.  To  get  a 
qualitative  look  at  this  effect,  analysis  with  the  electron  localization  function  (ELF)  is 
conducted. 
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Figure  29.  Logarithmic  plot  of  normalized  monomer/tetramer  ratio  for  seven 

different  Cp  variants. 
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d.  ELF  Analysis 

First  introduced  by  Becke  and  Edgecombe,  ELF  is  a  topological  analysis  of  a 
molecule  based  on  likely  locations  of  electron  pairs  [101].  ELF  takes  advantage  of  the 
Pauli  exclusion  principle  which  prevents  electrons  of  the  same  spin  from  occupying  the 
same  region  of  space  [102].  Using  the  electron  pair  density  function,  which  determine  the 
probability  of  finding  two  electrons  of  the  same  spin  in  given  regions  of  space,  one  can 
determine  where  pairs  of  electrons  of  opposite  spins  are  localized.  Because  of  the  Pauli 

exclusion  principle,  a  high  probability  of  finding  two  electrons  of  the  same  spin 

62 


corresponds  to  regions  where  electrons  are  delocalized  [103].  Conversely,  a  low 
probability  of  finding  two  electrons  of  the  same  spin  corresponds  to  regions  where 
electrons  are  well  localized  [103].  For  a  closed  shell  system  where  all  molecular  orbitals 
are  doubly  occupied,  as  is  the  case  for  all  clusters  studied  in  this  work,  regions  of  high 
electron  localization  would  correlate  to  locations  of  electron  pairs  of  opposite  spins.  ELF 
results  are  scaled  such  that  they  run  from  zero  to  one.  High  values  correlate  to  regions  of 
high  electron  pair  localization  and  low  values  corresponds  to  regions  of  low  pair 
localization.  A  value  of  0.5  corresponds  to  the  localization  of  electrons  in  a  uniform 
electron  gas.  The  results  help  to  connect  the  electron  density  with  Lewis’  concept  of  the 
electron  pair.  Figure  30  shows  the  ELF  results  for  the  AlCp  monomer  with  a  0.75 
isosurface  displayed.  The  graphic  gives  an  intuitive  picture  of  the  behavior  of  the  electron 
pairs  and  their  involvement  with  chemical  bonding.  The  red,  blue  and  green  surfaces 
display  behavior  of  valence  electrons  and  the  cyan  surfaces  show  core  shells.  The  red 
surface  corresponds  to  er  bonds  connecting  hydrogen  and  carbon  atoms.  The  green 
surface  corresponds  primarily  to  n  bonds  between  carbons,  but  also  includes  a  single 
electron  from  the  aluminum  atom.  The  blue  surface  corresponds  to  the  lone  pair  from  the 
aluminum  atom.  Finally,  the  cyan  surfaces  close  to  the  atomic  centers  are  core  shells.  The 
carbon  atoms  have  a  single  core  shell,  the  K  shell.  The  aluminum  atom  has  two  core 
shells,  the  K  and  L  shells. 
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Figure  30.  ELF  for  AlCp  monomer  with  the  0.75  isosurface  displayed. 

The  ELF  for  the  Cp-type  tetramers  is  very  similar  for  the  ligand  groups  with  clear 
bonding  regions  around  hydrogen  atoms  and  between  carbon  atoms.  The  ELF  in  the  core 
is  more  complex.  Figure  31  shows  the  0.79  isosurface  ELF  data  for  the  Cp,  CpA,  Cp*, 
Cp  Pr,  Cp  Ph  and  Cp  lPr  variants  spherically  clipped  so  that  only  data  within  4.5  A  of  the 
center  of  the  aluminum  tetrahedron  core  is  visible.  The  ELF  surfaces  for  the  two 
aluminum  core  shells  have  also  been  excluded.  Included  are  average  Al-Al  bond 
distances  for  each  variant.  For  a  better  understanding  of  the  significance  of  the  surface 
shapes,  ELF  bifurcation  analysis  is  conducted. 
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Figure  31.  ELF  for  aluminum  core  of  six  Cp-type  variants  with 
the  0.79  isosurface  displayed. 


In  order  to  understand  bifurcation,  it  is  necessary  to  understand  the  concept  of  an 
attractor  and  a  basin  in  ELF.  Attractors  are  essentially  maxima  on  the  ELF  surface. 
Basins  are  related  to  attractors  and  are  defined  as  the  region  of  space  within  which  the 
gradient  of  the  ELF  surface  points  towards  a  given  attractor.  In  Figure  30  there  are  17 
attractors  and  17  basins  surrounding  those  attractors.  Bifurcation  is  the  study  of  so-called 
/-localization  domains,  which  are  volumes  enclosed  by  ELF  isosurfaces  of  a  given  value, 
y  [104].  Figure  32(a)  depicts  the  ELF  of  C2H4  with  an / value  of  0.05  which  is  relatively 
close  to  zero.  At  this  level,  there  is  only  one  domain.  It  encloses  all  eight  attractors  in  the 
molecule:  the  two  carbon  core  attractors  and  the  six  valence  attractors  corresponding  to 

the  four  o  bonds  and  the  two  attractors  associated  with  the  double  n  bond.  As  the  / value 
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is  increased,  the  domain  will  shrink  and  will  eventually  bifurcate.  In  the  case  of  C2H4  the 
first  bifurcation  is  at  about  /=  0.15  at  which  point  the  single  domain  splits  into  three 
domains,  Figure  32(b).  Two  domains  are  considered  irreducible  because  they  enclose 
only  one  attractor,  the  center  of  each  carbon  atom.  The  third  domain  is  considered 
reducible  because  it  encloses  more  than  one  attractor  [105].  The  next  bifurcation  is  when 
/=  0.66  and  corresponds  to  the  point  when  the  remaining  reducible  domain  splits  into 
five  domains,  Figure  32(c).  Four  domains  contain  the  four  H-C  valence  attractors  and  the 
fifth  domain  contains  the  two  attractors  associated  with  the  C-C  double  bond.  The  final 
bifurcation,  not  shown  in  Figure  32,  occurs  at/=  0.93  when  the  domain  containing  the 
two  attractors  associated  with  the  C-C  double  bond  splits  in  two.  Bifurcation  analysis 
helps  understand  the  relative  electron  localization  between  basins.  The  lower  the 
bifurcation  point  for  a  given  attractor,  the  more  localized  the  electrons  are  within  the 
associated  basin  [104]. 


(a)  (b)  (c) 

Isosurface  values:  (a)/=0.05;  (b) /=  0.15;  (c)/=0.66. 

Figure  32.  Three  ELF  isosurface  values  illustrating  bifurcation. 


Figure  31  shows  aluminum  core  bonding  at  various  bifurcation  levels.  The  core  of 
the  Cp  variant  shows  six  clearly  distinguishable  domains  enclosing  the  six  Al-Al  valence 
attractors.  These  domains  are  all  irreducible  and  indicate  a  high  level  of  localization 
compared  with  the  other  five  variants.  This  would  suggest  that  the  Cp  variant  has  the 
strongest  Al-Al  core  bonding.  The  Cp  compound  shows  a  slight  merging  of  two  domains 
creating  a  single  reducible  domain  that  indicates  a  slight  electron  pair  delocalization  as 
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compared  with  the  Cp  variant.  This  is  interesting  considering  that  the  Cp  and  Cp  variants 
have  virtually  the  same  average  Al-Al  distance.  Insight  gained  through  ELF  could 
suggest  that  the  Al-Al  bonding  in  the  core  of  the  Cp  variant  is  marginally  weaker  than 
that  of  the  Cp  variant,  despite  similar  Al-Al  distances.  This  is  consistent  with  QTAIM 
analysis  which  shows  slightly  more  pb  associated  with  Al-Al  BCPs  in  Cp  variant.  The 
Cp  and  Cp*Pr  variants  show  additional  signs  of  delocalization  as  compared  to  the  Cp  and 
Cp  variants.  Both  have  similar  average  Al-Al  distances  and  both  show  the  presence  of 
reducible  domains.  There  are  two  reducible  domains  present  in  the  core  of  the  Cp*  variant 
and  only  one  in  the  core  of  the  Cp  Pr  variant.  The  two  domains  on  the  right  of  the  core  of 
the  Cp  Pr  variant,  however,  are  on  the  verge  of  merging  and  will  in  fact  merge  if  /  is 
lowered  to  0.78.  Additionally,  the  irreducible  domains  show  signs  of  distortion  as 
compared  to  the  Cp  variant,  indicating  some  electron  pair  delocalization.  The  Cp*Ph  and 
Cp  lPr  variants  have  the  largest  Al-Al  distances  of  any  of  the  variants  and  also  show  the 
most  amount  of  distortion  of  Al-Al  valence  domains.  The  Cp*Ph  variant  has  a  reducible 
domain  containing  three  attractors  and  though  the  three  irreducible  domains  at  the  bottom 
right  of  the  core  have  not  merged  yet,  they  will  at  /  =  0.78.  The  Cp*lPr  shows  two 
reducible  and  two  irreducible  domains,  similarly  to  the  Cp*  variant,  but  with  significant 
distortion  of  all  domains  as  compared  to  the  Cp*  variant.  These  two  variants  clearly  have 
the  most  amount  of  electron  delocalization  of  the  six  clusters  examined  and  therefore 
most  likely  have  the  weakest  Al-Al  bonding  in  the  core.  This  ELF  analysis,  while 
qualitative,  provides  insight  into  metal  bonding  in  the  core  that  compliments  QTAIM. 

Analysis  of  domains,  particularly  at  low  /  values,  has  also  been  used  to  study 
phenomenon  not  directly  captured  by  ELF  such  as  NCIs.  Covalent  bonds  between  atoms 
are  easily  identified  with  ELF  due  to  the  existence  of  an  attractor  between  atoms,  as  is  the 
case  for  C2H4.  The  enclosed  volume  between  the  two  carbon  atoms  in  Figure  32(c) 
encloses  the  valence  attractors  associated  with  the  double  n  bond  in  the  molecule.  NCIs, 
however,  cannot  be  identified  by  attractors  as  none  will  exist  between  the  interacting 
atoms.  Work  by  Massobrio  et  al.  showed  that  NCIs  could  be  identified  by  examining 
where  distinct  domains  join  as /is  decreased  [106].  They  showed  that  NCIs  between  OH' 
and  NO3'  in  copper  hydroxonitrate  could  be  identified  by  the  joining  of  the  domains  of 
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each  molecule  at /  values  of  about  0.15.  To  illustrate  this  point,  ELF  calculations  were 
conducted  on  a  typical  /?-tum  of  a  protein.  Proteins  are  long  strings  of  amino  acids,  but 
here  only  the  turn  is  investigated  and  connections  to  the  remainder  of  the  protein  chain 
were  truncated  with  hydrogen  atoms.  Figure  33  shows  the  ELF  of  the  /?-turn  at  an/ value 
of  0.06.  There  is  only  one  domain  encompassing  all  attractors  in  the  molecule,  but  the  left 
and  right  regions  of  the  domain  join  precisely  between  an  oxygen  atom  and  a  hydrogen 
atom.  This  identifies  a  known  NCI  that  is  the  primary  reason  the  /?-tum  exists  and 
ultimately  results  in  the  common  folding  behavior  of  large  proteins  [107]. 


Figure  33.  ELF  0.06  isosurface  showing  NCI  between  hydrogen  and  oxygen 

atoms  in  a  /?-turn  of  amino  acids. 

The  connecting  of  separate  domains  as  well  as  the  connecting  of  separate  regions 
of  the  same  domain  both  occur  within  all  Cp-type  variants.  Figure  34  shows  two 
neighboring  ligands  of  the  Cp*Ph  variant  with  most  of  the  ELF  isosurface  clipped  out  for 
clarity.  The  isosurface  corresponds  to /=  0.03.  Highlighted  with  the  red  circle  is  an  inter¬ 
ligand  connection  between  the  separate  domains  of  the  two  ligands.  This  connection  is 
associated  with  an  NCI  between  a  Me  hydrogen  atom  from  the  lower  ligand  and  a  Me 
carbon  atom  from  the  upper  ligand.  This  also  correlates  with  an  H-C  BCP  discovered 
through  QTAIM  analysis.  Highlighted  with  green  circles  are  intra-ligand  connections 
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between  two  regions  of  the  domains  surrounding  each  ligand.  The  connections  are 
between  Me  hydrogen  atoms  and  Ph  carbon  atoms  from  the  same  ligand.  Both  also 
correlate  with  QTAIM  H-C  BCPs.  This  would  suggest  that  connections  between  separate 
domains  or  between  separate  region  of  the  same  domain  can  indicate  bone-fide  NCIs. 


Figure  34.  ELF  0.03  isosurface  showing  inter-ligand  connection  (red  circle)  and 
intra-ligand  connections  (green  circles)  of  domains  in  Cp*Ph  variant. 

2.  Other  Isopropyl  Variants 

Because  the  addition  of  bulk,  in  general,  does  not  to  lead  to  a  higher  monomer  to 
tetramer  ratio  for  Cp-type  clusters,  three  variants  that  have  functional  groups  which  may 
hinder  bonding  in  the  aluminum  core  are  examined.  Investigation  of  these  variants  was 
suggested  by  researchers  at  the  China  Lake  NAWC.  They  are  based  on  the  Cp*lPr  variant 
and  all  have  the  form  AlCpniPl  (CpniPr  =  C5iPrnH(5.n);  n=3,4,5).  The  Cp*lPr  variant  has  the 
largest  Al-Al  average  bond  distance  in  the  core  of  the  tetramer  and  despite  being  one  of 
the  most  enthalpically  stable  variants,  it  has  the  highest  experimental  monomer  to 
tetramer  ratio  of  any  variant  previously  discussed.  All  three  CpmPr  variants  were  subjected 
to  the  same  theoretical  structural  and  thermochemical  analysis  as  the  Cp  Ph  and  their 
potential  monomer/tetramer  ratios  calculated. 
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a.  Triisopropyl  Variant 

Before  attempting  to  find  the  lowest  energy  tetrameric  configuration  for  the  Cp3lPr 
variant,  DFT  analysis  into  the  monomer  was  conducted.  There  are  two  stable  directions 
the  lone  hydrogen  atom  of  the  iPr  functional  group  can  point  in  the  monomer.  The 
monomer,  aluminum  (I)  1,2,4-triisoproylcyclopentatienide,  has  iPr  groups  on  the  ligand 
at  the  1,  2  and  4  positions.  For  the  iPr  group  at  position  4,  a  change  in  direction  is  has  no 
impact  on  ground  state  energy.  This  leaves  four  possible  combination  of  directions  for  the 
two  remaining  groups.  Due  to  symmetry  only  three  of  these  configurations  are  distinct. 
Ground  state  energies  for  all  three  types  were  calculated  with  M06-2X/cc-pVTZ  and  the 
configuration  shown  in  Figure  35  was  found  to  be  6.5  kJ/mol  lower  in  energy  than  the 
next  closest  configuration.  For  this  reason,  tetramers  created  with  this  monomer  were 
ultimately  used  for  the  GA. 


Figure  35.  Lowest  energy  configuration  of  the  AlCp3lPl  monomer. 

The  lowest  energy  structure  for  the  Cp3lPr  tetramer  calculated  via  the  GA  found 
the  lowest  structure  in  the  9th  generation,  though  this  structure  was  only  0.13  kJ/mol 
lower  in  energy  than  the  previous  lowest  structure  discovered  in  the  5th  generation.  The 
minimum  found  by  the  GA  was  relaxed  using  M06-2X/cc-pVTZ  in  Gaussian. 
Surprisingly,  the  average  Al-Al  distance  for  this  variant  was  2.746  A,  making  it  slightly 
more  compact  than  the  Cp  variant  and  significantly  more  compact  than  the  Cp  lPr 


70 


variant.  The  inter-ligand  dispersive  forces  due  to  the  extra  bulk  are  likely  strong  enough 
to  overcome  the  hindrance  of  the  iPr  functional  groups  and  pack  the  aluminum  core 
closer  than  expected.  This  produces  a  synergistic  effect  whereby  the  additional  stability 
gained  through  inter-ligand  NCIs  actually  forces  aluminum  atoms  closer  which,  in  turn, 
increases  the  stability  gained  through  Al-Al  bonding.  This  is  reflected  in  the  calculated 
A Htet  value  of  -181.1  kJ/mol.  This  makes  it  slightly  less  stable  then  the  Cp*Ph  variant 
which  has  an  A Htet  value  of  -182.9  kJ/mol. 

ELF  analysis  of  the  Cp3lPr  variant  indicate  very  localized  bonding  in  the  aluminum 
core.  Figure  36  shows  the  Al-Al  valence  ELF  isosurface  with/=  0.79.  The  Cp3lPr  variant 
has  six  distinct  irreducible  domains  surrounding  the  six  Al-Al  bond  attractors.  The  core 
of  this  variant  looks  very  similar  to  the  core  of  the  Cp  variant  and  more  localized  than 
even  the  Cp  which  has  a  smaller  average  Al-Al  distance.  It  would  appear  that  despite 
having  a  larger  average  Al-Al  distance,  the  aluminum  bonding  in  the  core  of  the  Cp3lPr 
variant  may  be  stronger  than  that  of  the  Cp  variant. 


Figure  36.  ELF  for  aluminum  core  of  the  Cp3lPr  tetramer  with 
the  0.79  isosurface  displayed 
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b.  Tetraisopropyl  Variant 

There  are  several  different  configurations  of  the  Cp4lPr  monomer.  The  four  iPr 
groups  can  be  arranged  in  16  different  configurations.  Due  to  symmetry,  only  eight  are 
distinct.  The  ground  state  energies  for  all  eight  configurations  were  calculated  with 
M062X/cc-pVTZ.  The  lowest  energy  configuration  is  shown  in  Figure  37  and  was 
ultimately  used  to  create  tetramers  for  the  GA.  The  Cp4lPl  proved  to  be  much  more 
difficult  to  pack  together  then  the  Cp  ,lPr  variant.  After  ten  generations  of  the  GA,  the 
lowest  energy  structure  was  optimized  in  Gaussian  and  found  to  have  an  average  Al-Al 
bond  distance  of  3.4567  A,  making  the  Al-Al  core  bonding  significantly  weaker  than  all 
previous  variants. 


Figure  37.  Lowest  energy  configuration  of  the  AlCp4lPl  monomer. 

The  calculated  A Htet  for  the  Cp4lPr  variant  is  -72.9  kJ/mol,  making  this  by  far  the 
least  enthalpically  stable  variant.  This  result  is  primarily  due  to  the  increased  average  Al- 
Al  bond  distance.  QTAIM  analysis  indicates  that  the  Al-Al  bonds  here  are  distinct  from 
those  of  the  other  clusters  analyzed.  All  six  Al-Al  BCPs  are  associated  with  a  V2pb  >  0, 
which  is  unique  for  this  variant,  and  a  relatively  low  total  value  of  pb,  0.0840.  All  other 
variants  have  V2pb  <  0  and  an  average  total  pb  of  0.2127  at  Al-Al  BCPs.  A  V2pb  >  0 
combined  with  a  relatively  low  value  of  pb  at  a  BCP  is  an  indication  of  local  charge 
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depletion  which  is  typically  associated  with  NCIs  [108].  As  Bone  and  Bader  point  out, 
however,  a  V2p&  <  0  is  not  a  necessary  indicator  of  a  covalent  bond  [85].  The  Al-Al 
BCPs  for  the  Cp4lPr  variant  are  also  associated  with  Hb  <  0,  a  kinetic  energy  per  density 
ratio  ( Gb/Pb )  much  less  than  1  and  \Vh\  >  Gb.  These  are  characteristics  typically 
associated  with  covalent  bonding  [82].  The  unique  nature  of  this  bonding  is  also  reflected 
in  ELF  analysis  of  the  aluminum  core.  Figure  38  shows  core  ELF  isosurfaces  for  three 
different/ values.  With/=  0.76  it  is  clear  that  the  six  Al-Al  valence  attractors  present  in 
all  other  variants  are  non-existent.  The  surfaces  seen  in  Figure  38  are  those  enclosing  the 
four  attractors  associated  with  the  lone  pair  of  electrons  associated  with  each  aluminum 
atom.  Connections  between  the  four  domains  begins  as/is  reduced  through  0.47.  This  is 
seen  in  the  middle  image  of  Figure  38  where  the  domains  of  the  two  rearward  A1  atoms 
have  joined.  Connections  between  all  aluminum  atoms  is  achieved  when/=  0.36.  Each 
aluminum  domain  is  connected  to  the  domains  of  other  three  aluminum  atoms  creating 
six  distinct  connections  analogous  to  the  six  valence  attractors  present  in  each  of  the 
previous  eight  clusters  examined.  Connections  between  aluminum  domains  occurring 
between  /=  0.36  and  /=  0.47  are  below  the  typical  value  of  /  >  0.5  associated  with 
covalent  bonds  [109].  This,  in  addition  to  QTAIM  analysis,  suggests  that  describing 
aluminum  lone  pair  overlap  in  the  core  of  the  Cp4lPr  variant  as  covalent  is  tenuous. 


/=  0.76 


/=  0.47 


/=  0.36 


Figure  38.  ELF  for  aluminum  core  of  the  Cp4lPr  tetramer. 
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c.  Pentaisopropyl  Variant 

The  Cp5lPl  monomer  has  an  iPr  group  attached  to  every  carbon  site  on  the  Cp  ring. 
For  this  reason,  all  iPr  groups  must  be  orientated  in  the  same  direction.  Switching  the 
orientation  of  one  group  results  in  about  a  25  kJ/mol  increase  in  ground  state  energy,  and 
switching  multiple  groups  results  in  even  higher  energy  changes.  This  leaves  two 
possible  configurations  for  the  Cp"'1’1  monomer,  but  they  are  degenerate.  The 
configuration  used  for  to  build  tetramers  for  the  GA  is  shown  in  Figure  39.  Despite  a  run 
that  included  15  generations  and  more  than  200  structures,  the  GA  never  produced  a 
Cp5lPr  variant  with  an  average  Al-Al  distance  smaller  than  4.9  A.  This  is  well  outside  the 
Al-Al  bonding  region  and  all  subsequent  attempts  to  optimize  top  GA  candidates  in 
Gaussian  failed.  The  addition  of  the  fifth  iPr  group  to  each  monomer  appears  to  have 
provided  enough  hindrance  to  prevent  aluminum  lone  pair  overlap  and  as  a  consequence 
no  stable  tetramer  was  found.  While  synthesis  of  the  monomer  is  possible,  as  Fischer’s 
reductive  elimination  method  produces  monomers  directly  [63],  it  is  unclear  if  larger 
metalloid  clusters  would  be  stable  without  a  stable  tetramer  for  this  variant.  Much  is 
unknown  about  the  growth  process  of  clusters  in  solution,  however,  and  it  is  possible  that 
cluster  growth  routes  which  do  not  involve  the  tetramer  may  exist. 
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d.  Equilibrium  Constants  for  CpmPr  Variants 

The  three  iPr  variants  just  discussed  have  a  steady  progression  of  increasing  steric 
bulk  and  hindrance  preventing  short  Al-Al  bond  distance.  The  variant  with  the  smallest 
amount  of  bulk,  Cp3lPr,  is  the  most  enthalpically  stable  due  primarily  to  the  fact  that, 
despite  its  bulk,  it  has  an  average  Al-Al  bond  distance  even  smaller  than  those  of  the  Cp 
variants.  The  hindrance  of  Cp4lPr  and  Cp5lPr  variants,  however,  force  much  larger  Al-Al 
bond  distances  and  in  the  case  of  the  Cp5lPr  variant,  likely  prevent  the  existence  of  a 
stable  tetramer.  Figure  40  shows  a  logarithmic  plot  of  monomer/tetramer  ratios  for  the 
Cp,  Cp*,  Cp*Ph,  Cp3lPr  and  Cp4lPr  variants  normalized  so  that  the  ratio  for  the  Cp*  variant 
is  1:449,  the  known  experimental  ratio.  Results  indicate  that  the  monomer  to  tetramer 
ratio  for  the  Cp3lPr  variant  would  fall  between  that  of  the  Cp*  and  Cp  Ph  variants.  Whether 
or  not  monomers  for  this  variant  would  be  found  in  solution  is  inconclusive.  The  Cp4lPl 
variant,  however,  has  a  much  lower  calculated  equilibrium  constant  and  therefore  a  much 
higher  monomer  to  tetramer  ratio.  The  ratio  of  the  Cp4lPr  variant  is  actually  higher  than  all 
other  variants  examined  and  is  that  only  variant  with  more  monomers  than  tetramers  in 
solution.  This  likely  makes  it  the  best  precursor  for  larger  metalloid  clusters. 
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Figure  40.  Logarithmic  plot  of  monomer/tetramer  equilibrium  ratios  for  Cp,  Cp*, 

Cp*Ph,  Cp1lPr  and  Cp4lPr  variants. 
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IV.  CONCLUSION 


Large  metalloid  clusters  such  as  Al50Cp  12  show  enormous  potential  as  novel 
energetic  materials.  Theoretical  heats  of  combustion  suggest  that  these  clusters  have  very 
high  energy  density  compared  to  explosives,  but  their  small  size  and  core/shell  structure 
may  allow  for  much  more  rapid  combustion  kinetics  than  standard  metal  fuels.  In  this 
work,  computational  analysis  of  several  small,  prototypical  monovalent  aluminum 
clusters  was  perfonned  with  the  ultimate  goal  of  developing  novel  energetic  materials. 
DFT  calculations  show  that  the  predicted  thermochemistry  is  extremely  sensitive  to  the 
functional  choice,  both  due  to  the  prominence  of  non-covalent  interactions  between 
ligands  and  DFT  inaccuracies  in  treating  the  low-valent  A1  bonding  in  the  core. 
Calculated  results  compared  with  available  experimental  thermochemical  data  for  three 
different  AlCp-based  variants  indicate  the  necessity  to  account  for  short  and  medium 
range  dispersive  forces.  This  work  demonstrates  that  the  Minnesota  functional  M06-2X, 
specifically  parameterized  to  capture  these  effects,  is  an  exceptional  method  for 
calculating  key  characteristics  of  these  clusters.  It  predicted  enthalpies  of  tetramerization 
within  2%  of  experimental  values  and  geometrical  parameters  to  better  than  1%  of 
available  data  when  paired  with  the  triple  zeta,  cc-pVTZ  basis  set.  This  work  also 
demonstrates  that  M06-2X  is  more  accurate  and  computationally  expedient  than  methods 
that  used  additional  calculations  to  explicitly  treat  dispersive  interactions. 

Using  a  GA  and  M06-2X,  the  lowest-energy  structures  of  ten  AlCp-based  clusters 
with  varying  amounts  of  steric  bulk  were  determined.  Thermochemical  analysis  indicated 
the  addition  of  steric  bulk  on  the  Cp  ligands  tends  to  stabilize  these  clusters  enthalpically, 
primarily  via  NCIs  across  ligand  groups.  QTAIM  analysis  shows  that  there  is  a  strong 
correlation  between  inter-ligand  BCPs  and  overall  stability  of  the  tetrameric  form. 
Additionally,  for  all  variants  except  the  unsubstituted  AlCp  variant,  covalent  Al-Al 
bonding  in  the  tetrahedral  core  accounted  for  only  about  65%  of  the  overall  enthalpic 
stability  of  the  tetramer.  Equilibrium  analysis,  however,  reveals  no  trend  between  ligand 
steric  bulk  and  monomer  to  tetramer  equilibrium.  The  two  variants  with  the  least  amount 
of  steric  bulk  have  by  far  the  lowest  monomer  to  tetramer  ratio.  The  Cp  Ph  variant, 
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however,  has  the  lowest  ratio  of  monomer  to  tetramers  amongst  all  the  Cp  variants, 
despite  having  the  most  steric  bulk.  Additional  steric  bulk  seems  to  only  increase 
monomer  to  tetramer  ratio  when  the  bulk  interferes  with  Al-Al  covalent  bonding  in  the 
tetrahedral  core.  The  addition  of  a  single  iPr  group  to  the  Cp3lPr  variant  is  sufficient  to 
force  the  average  Al-Al  distance  further  apart  by  approximately  0.6  A.  The  resulting 
variant,  Cp4lPr,  has  a  tetramer  with  the  most  steric  bulk  of  any  examined,  but  also  with  the 
second  weakest  enthalpy  of  tetramerization.  This  reduction  in  stability  produces  the 
highest  monomer  to  tetramer  ratio  of  all  variants,  making  it  a  good  candidate  for  solution- 
based  synthesis  techniques. 

DFT  with  M06-2X  revealed  that  there  is  no  trend  relating  bulk  and  monomer  to 
tetramer  ratio  in  solution.  Each  potential  cluster  must  be  examined  individually.  The  hope 
is  that  cluster  variants  with  large  numbers  of  monomers  in  solution  will  provide  the 
building  blocks  needed  to  produce  larger  systems  analogous  to  Al50Cp*i2.  Other  avenues 
to  large  metalloid  clusters  are  likely  possible  and  should  be  pursued.  The  ultimate  goal  is 
to  passivate  aluminum  atoms  with  specifically  designed  organic  ligands  to  create  clusters 
with  accelerated  combustion  kinematics  compared  to  much  larger  aluminum  particles  that 
are  typically  added  to  high  explosives. 

Beyond  metalloid  cluster  synthesis,  it  is  also  import  to  continue  to  examine 
potential  stabilizing  ligands.  The  ligands  examined  in  the  work  contain  no  oxidizing 
agents.  Aluminum  combustion  would  require  the  addition  of  an  oxidizer  to  the  eventual 
explosive  compound  produced.  Ligands  whose  bulk  includes  multiple  nitro  groups,  for 
example,  would  allow  for  the  inclusion  of  some  oxidizing  agent  as  part  of  the  cluster. 
This  would  reduce  the  amount  of  external  oxidizer  required.  Synthesizing  metalloid 
clusters  passivated  with  fluorinated  monomers  such  as  AlCp  5F  (Cp  31  =  C5(CF3)5)  could 
create  material  with  enough  oxidizer  for  the  entire  core  of  aluminum  atoms,  potentially 
eliminating  the  need  for  an  added  oxidizer.  Theoretical  analysis  of  these  and  other  types 
of  ligands  will  help  determine  their  effect  on  tetrameric  stability  as  well  as  the  viability  of 
these  tetramers  to  produce  larger  clusters. 


78 


APPENDIX  A.  COMPUTATIONAL  METHODOLOGY 


DFT  calculations  of  all  compounds  were  carried  out  using  the  Gaussian  G09.E01 
and  G16.A03  software  [110],  [111].  All  calculations  used  an  ultra-fine  grid  and  systems 
were  relaxed  until  no  imaginary  vibrational  frequencies  were  found.  Basis  set 
superposition  error  analysis  was  done  using  a  counterpoise  methodology  [112].  Solvent 
corrections  were  made  using  the  conductor-like  screen  model  (COSMO)  as  implemented 
in  the  Polarizable  Continuum  Model  (PCM)  framework  for  Gaussian  [113],  [114]. 
Calculations  which  include  explicit  dispersion  corrections  were  done  using  Grimme’s  D3 
approach  [54].  Bond  critical  points  and  associated  density  metrics  were  determined  using 
the  DGrid  software  [115].  In  order  to  calculate  27A1  NMR  chemical  shifts,  known 
experimental  NMR  values  were  plotted  against  the  isotropic  output  from  a  GIAO 
calculation  scaled  by  the  0  ppm  reference,  A1(H20)63'  [116].  ELF  analysis  was  done 
using  the  TopMoD  software  [1 17] — [1 19].  All  calculations  were  performed  on  DOD  HPC 
computer  clusters. 
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APPENDIX  B.  ADDITIONAL  PUBLICATIONS 


The  following  publications  resulted  from  this  dissertation  work: 


S.  Vaddypally,  I.  G.  McKendry,  W.  Tomlinson,  J.  P.  Hooper,  and  M.  J.  Zdilla, 
“Electronic  structure  of  manganese  complexes  of  the  redox-non-innocent 
tetrazene  ligand  and  evidence  for  metal-azide/imido  cycloaddition  intermediate,” 
Chem.  Eur.  J.,  vol.  22,  pp.  10548-10557,  2016. 

S.  DeCarlo,  D.  H.  Mayo,  W.  Tomlinson,  J.  Hu,  J.  P.  Hooper,  P.  Zavalij,  K. 
Bowen,  H.  Schnockel,  and  B.  Eichhom,  “Synthesis,  structure  and  properties  of 
Al(Rbpy)3  complexes  (R=/-Bu,  Me):  Homoleptic  maingroup  tris-bipyridyl 
compounds,”  Inorg.  Chem.,  vol.  55,  pp.  4344-4353,  2016. 

J.  R.  Heinmer,  S.  O.  Peolrna,  N.  Treat,  Z.  A.  Page,  N.  D.  Dolinski,  Y.  J.  Diaz,  W. 
Tomlinson,  K.  D.  Clark,  J.  P.  Hooper,  C.  Hawker,  and  J.  R.  Alaniz,  “Tunable 
visible  and  near  infrared  photoswitches”,  J.  Am.  Chem.  Soc.,  vol.  138,  pp.  13960- 
13966,2016. 

W.  Tomlinson,  D.  H.  Mayo,  R.  M.  Wilson,  and  J.  P.  Hooper,  “The  role  of  steric 
bulk  in  new  monovalent  aluminum  compounds,”  J.  Chem.  Phys.  A.,  vol.  121,  no. 
24,  pp.  4678-4687,2017. 
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